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spots,  is  with  anode  gas  injection  clearly  distinct  from  it.  The 
incidence  of  these  instabilities  is  also  delayed  with  rising  an¬ 
ode  gas,  but  in  a  slower  manner  than  the  spot  phenomenon.  This 
indicates  different  onset  mechanisms.  Nitrogen  as  propellant 
shows  no  merit  over  argon  in  the  investigated  current  ranges; 
the  discharge  voltage  is  higher  and  thus  it  lowers  the  efficien¬ 
cy.  First  experiments  with  a  cylindrical  thruster  are  demon¬ 
strated.  They  show  the  feasibility  of  a  hot,  protruding  cathode 
and  anode  current  attachment  dependencies  on  the  cathode  posi¬ 
tion.  A  concept  of  a  200  -  300  kW  MPD  thruster  with  a  radiation- 
cooled  anode  is  presented. 

Since  there  are  differences  in  the  results  of  steady  state  and 
quasi-steady  pulsed  thrusters,  experiments  with  thrusters  simi¬ 
lar  in  geometry  and  electrode  material  have  also  been  performed 
in  quasi-steady  operation  mode.  They  indicate  lower  performances 
due  to  lower  thrust,  higher  discharge  voltages  and  lower  onset 
points.  Data  o£  quasi-steady  thrusters  cannot  be  extrapolated  to 
the  continuous  mode. 

Cathode  erosion  studies  were  made  for  cold  cathodes  (pulsed 
thrusters  and  ignition  phase  of  steady  state  ones)  and  hot  cath¬ 
odes.  They  indicate  an  erosion  of  almost  three  orders  of  magni¬ 
tude  higher  with  the  cold  cathode  than  with  the  hot  ones.  An  ex¬ 
planation  by  spotty  arc  attachment  to  the  cathode  is  given. 

Three  theoretical  models  for  calculating  the  current  and  flow 
distributions  in  self-field  MPD  thrusters  are  being  developed 
and  will  be  applied  to  evaluate  the  effects  of  geometry,  propel¬ 
lant  type,  scale  and  other  parameters  on  the  thruster  perform¬ 
ance  . 

1)  A  sem  i- anal yt ica 1 ,  two-dimensional  code  based  on  a  source 
flow  model. 

2)  For  continuous  thrusters,  a  stationary  code  is  being  devel¬ 
oped.  The  extended  Ohm's  law  is  used  to  calculate  the  current 
contour  lines,  and  a  one-dimensional,  two-component  expansion 
flow  model  is  used  to  obtain  the  velocity,  temperature  and 
pressure  distributions  for  calculating  the  gas  properties, 
which  are  again  used  in  Ohm’s  law.  The  differential  equation 
is  solved  by  means  of  a  finite  difference  method  for  the  geo¬ 
metry  of  the  nozzle-type  plasma  thruster  DT2 ,  which  has  been 
investigated  at  LkS  in  a  steady  state  as  well  as  in  a  quasi- 
steady  state  mode. 

An  integration  over  the  volume  and  thermal  forces  equals  to 
the  thrust.  The  calculated  current  density  distribution  and 
tne  computed  thrust  are  compared  with  experimental  results. 

J)  For  the  starting  phase  of  the  steady  state  MPD  thrusters  as 
well  as  for.  pulsed  tnrusters,  a  time  dependent,  fully  two- 
dimensional  code  is  being  developed.  it  uses  a  modified 
MeCo  t  nine  K  FD  method  in  cylindrical  coordinates  to  calculate 
the  time  dependent  flow,  temperature  and  pressure  fields.  The 
current  and  magnetic  contour  lines,  which  must  be  known  for 
tne  energy  and  impulse  input  in  this  code,  are  calculated 
witn  an  instat  ionary  Ohm's  law  solved  in  a  manner  similar  to 
the  computer-  cntP  for  continuous  thrusters.  This  code  is  useu 
for  .1  simplified  cyl  nnlru’il  thruster  'node  1  . 
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1 .  Introduction 

Despite  the  present  great  lead  of  ion  engines  in  development 
level,  performance  and  life  for  sizes  up  to  several  kW,  new  in¬ 
terest  in  magnetoplasmadynamic  thrusters  for  higher  power  levels 
(up  to  hundreds  of  kW)  has  arisen  in  recent  years.  Mission  stud¬ 
ies1,2  showed  specific  impulse  values  between  1000  and  2000  s 
for  orbit  raising  and  orbit  transfers  for  optimized  mass  savings 
and  trip  times,  and  advantages  for  MPD  propulsion  even  with  rel¬ 
atively  low  efficiencies. 

With  the  probable  realization  of  the  SP-100  nuclear  space  power 
plant,  sufficient  power  would  also  be  available  for  steady  state 
MPD  propulsion,  so  that  the  benefits  of  these  thrusters  versus 
pulsed  ones  could  be  utilized;  i.e.,  the  ruggedness  and  simplic¬ 
ity  of  power  and  propellant  conditioning  and  the  about  3  orders 
of  magnitude  smaller  cathode  erosion.  But  there  are  also  disad¬ 
vantages  of  steady  state  thrusters;  e  .g . ,  with  quasi-steady  ones 
(but  with  other  geometries)  the  achieved  efficiencies  and  spe¬ 
cific  impulses  were  in  some  cases  higher.  J 

The  research  work  covered  in  this  report  deals  with  these  and 
related  problems:  in  chapter  two,  the  experimental  results  with 
nozzle- type  and  cylindrical  thrusters  are  presented.  The  influ¬ 
ence  of  cathode  position,  ambient  pressure,  anode  gas  injection 
and  propellant  type  on  the  performance  is  investigated.  The 
thruster  design  studies  and  the  heat  transfer  calculations  for 
the  MPD  thruster  with  a  radiation-cooled  anode  are  reported. 

With  the  continuously  running  Stuttgart  thrusters,  the  critical 
iVm  values  and  therefore  the  attainable  specific  impulses  are 
smaller  compared  with  the  data  read  at  Princeton  and  Tokyo.  To 
check  whether  these  differences  are  a  function  of  the  design 
configuration  or  of  the  operation  mode,  thrusters  identical  in 
geometry  and  electrode  material  were  built  and  tested  in  both 
the  continuously  running  and  the  pulsed  mode.  The  experimental 
results  are  discussed  in  chapter  3. 
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Cnapter  4  reports  about  erosion  measurements  on  pulsed  and  con¬ 
tinuous  MPD  thrusters,  in  the  case  of  steady  state  operation, 
two  different  kinds  of  erosion  must  be  distinguished:  erosion 
during  starting  phase,  when  the  cathode  is  cold,  and  the  much 
lower  erosion  during  steady  state  operation,  when  the  cathode  is 
hot.  For  both  phases,  the  cathode  mass  loss  has  been  separately 
measured.  Erosion  measurements  on  a  pulsed  thruster  with  the 
same  electrode  material  and  geometry  allow  the  comparison  of 
both  operationmodes  with  respect  to  electrode  erosion. 

The  theoretical  evaluation  of  the  electrode  erosion  mechanisms, 
also  a  part  of  chapter  4,  has  been  carried  out  under  contract 
AFOSR  82-0298. 4 

Parallel  to  the  experimental  investigations,  analytical  and  nu¬ 
merical  approaches  to  determine  flow,  pressure,  temperature  and 
current  density  field  of  the  midstream  flow  discharge  region  in 
MPD  thrusters  were  carried  out  (partly  under  Grant  AFOSR  82- 
0298).  These  programs  are  applied  to  the  specific  geometries  of 
the  Stuttgart  thrusters.  In  chapter  5,  results  are  demonstrated 
and  compared  with  experiments. 

The  steady  state  experiments  presented  are  made  possible  by  a 
facility  unique  in  Europe  with  highly  regulated  (low  ripple)  dc 
power  up  to  48,000  A  and  6  MW  continuous.  A  large  flexible  vacu¬ 
um  system  with  up  to  about  33  m3/s  (33,000  1/s)  pumping  capacity 
permits  continuous  testing  over  a  wide  range  of  tank  pressures 
and  massflow  values.  Tank  pressures  at  or  below  2  •  1 0 —  mbar , 
adequate  at  least  for  self-field  thrusters,  can  be  maintained 
with  a  flow  rate  of  up  to  1  g/s  with  argon.  The  vacuum  system 

__  I4 

olank-off  pressure  (without  the  large  tanks)  is  4  x  10  mbar, 
and  with  tanks  connected  at  present  about  10”  mbar--the  dif¬ 
ference  being  due  to  moisture  and  other  outgassing,  e  .g .  from 
leaks  and  from  the  silicone  oil  of  the  thrust  balance. 

In  Fig.  1  a  view  of  the  laboratory  hall  of  the  IRS,  as  of  summer 
1986,  is  shown  with  the  different  vacuum  tanks  and  the  computer 
plat  form. 


2.  Steady  State  MPD  Thrusters 


2. 1  Nozzle-Type  Thruster 

The  emphasis  of  the  design  and  test  work  during  the  research 


period  has  been  on 

nozzle-type  self-field 

MPD 

thrusters . 

Two  different  design 

concepts  with  both  the 

same 

nozzle 

and 

electrode  dimensions 

were  investigated.  While 

the 

first , 

DTI 

(Fig.  2),  originates  from  a  DFVLR  concept  of  Hugel  5  and  was 
built  for  argon  as  propellant,  the  newer  DT2  (Fig.  3),  with 
thinner  segmentation  allowed  molecular  propellants  as  well  as 
argon.  Thinner  segmentation  avoids  breakdown  between  the  seg¬ 
ments  and  is  necessary,  since  the  voltage  drop  along  the  arc 
axis  and  the  nozzle  wall  increases  by  taking  nitrogen  and  hydro¬ 
gen  as  propellants. 

The  thrusters  consist  of  an  upstream  cathode  made  of  2%  thori- 
ated  tungsten  which  is  cast  in  a  water-cooled  copper  rod,  an  arc 
chamber,  and  an  expansion  nozzle.  The  arc  chamber  and  nozzle  are 
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made  of  individual  water-cooled  segments;  these  are  insulated 
against  each  other  and  are  electrically  floating.  The  last  and 
largest  of  the  segments  represents  the  anode. 

The  thruster's  dimensions  include  a  throat  diameter  of  24  mm 
and  an  outer  anode  diameter  of  142  mm.  All  parts  of  the  anode 
with  greater  diameter  are  covered  with  a  ceramic  insulation,  in 
this  manner  the  region  of  the  arc's  anode  attachment  is  rela¬ 
tively  well  defined.  The  propellant  is  fed  mainly  through  the 
nozzle  chamber,  but  additional  gas  can  be  injected  at  the  anode 
ring  with  the  DT 2  thruster. 

The  thrusters  are  operated  at  currents  of  up  to  6000  A  with  pro¬ 
pellant  flows  between  0.2  and  1.5  g/s,  mainly  using  argon. 

For  direct  thrust  measurements,  the  thruster  is  mounted  on  a 
thrust  balance  and  placed  into  a  vacuum  tank  with  a  length  of 
about  5  m  and  a  diameter  of  2  m  and  which  has  double,  water- 
cooled  walls.  The  current  feed  to  the  balance  is  achieved 
through  mercury  baths.  Propellant  and  cooling  water  are  supplied 


insulator 


Fig.  3  Nozzle-type  MPD _ thruster  DT2 _ with  anode  gas  injection 

e_<iu  ipinent . 
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by  nylon  hoses  which  form  large  bows  in  vertical  planes.  In 
order  to  eliminate  errors  due  to  electromagnetic  and  mechanical 
interactions,  calibration  measurements  have  been  performed  with 
the  thruster  short  circuited  and  with  running  cooling  water.  In 
addition,  the  thrust  balance  is  remotely  calibrated  always  be¬ 
fore  and  after  every  test,  allowing  thrust  measurements  within 
an  error  of  less  than  2%.  The  schematic  of  the  thrust  stand  is 
shown  in  Fig.  4.  Further,  the  current  discharge  voltage, 
propellant  mass  throughflow,  and  thermal  losses  for  each 
thruster  segment  and  for  the  electrodes  are  measured. 


The  measured  dates  were  recorded  by  an  automatic  data  acquisi¬ 
tion  system  which  is  connected  to  a  DEC  laboratory  computer. 


Experimental  results 


Influence  of  the  ambient  pressure  on  thruster  performance 
(measured  with  DTI  thruster) 

In  some  recent  papers,*3  questions  were  raised  about  the  validity 
of  test  results  of  continuous  self-field  MPD  thrusters  because 
of  the  relatively  high  ambient  pressure  of  about  5*10"J  to 
5*10"2  mbar  in  the  vacuum  systems.  These  authors  referred  to 
measurements  which  were  not  made  with  self-field,  but  with  axial 
appl ied- f ield  ( MPD-A  or  Hallarcjet  or  MAARC)  thrusters.  This 
sort  of  thruster  indeed  showed  a  strong  influence  by  the  ambient 
pressure  on  its  per formance :  ' , 8 , 9  the  curve  of  measured  thrust 
versus  vacuum  tank  pressure  showed  a  declining  thrust  with  in¬ 
creasing  pressure  in  the  region  between  to-4  and  10"  1  mbar. 

But  there  is  no  indication  anywhere  in  the  literature  that  any 
tests  of  continuous  self-field  thrusters  at  tank  pressures  in 
the  regime  reported  here  (~  1-2.5x10" 2  mbar)  have  produced  erro¬ 
neous  or  questionable  data.  This  has  been  confirmed  by  measure¬ 
ments  in  our  facility,  shown  in  Fig.  5.  With  increasing  tank 
pressure,  there  was  initially  no  measurable  change  in  thruster 
performance.  Only  large  changes  in  pressure  (by  a  factor  of  ~  5 


. 


Locking  and  Calibrating  Mechanism 


there  appears  to  be  no  possibility  that  tank  gas  either  gets 
sucked  in  or  diffuses  into  that  part  of  the  discharge  which  de¬ 
termines  the  measured  performance.  This  is  particularly  valid 
for  steady  state  thrusters,  which  are  generally  operated  at 
higher  than  the  critical  mass  flow,  or  at  I 2/m  well  below  onset 
value.  Far  out  in  the  plume,  the  tank  gas  of  course  mixes  with 
the  jet  and  retards  the  thruster  plasma  plume  (as  Kuriki 1 1  and 
others  have  measured),  but  this  has  no  influence  on  the  measured 
performance . 


Effect  of  cathode  position  (measured  with  DTI  thruster) 

To  determine  the  effect  of  the  cathode  position  on  the  thruster 
performance,  several  tests  with  varying  cathode  positions  were 
made,  in  the  course  of  which  a  strong  effect  on  the  thrust,  the 
arc  voltage  and  therefore  also  on  the  efficiency  had  been  no¬ 
ticed.  With  the  cathode  tip  closer  to  the  nozzle  throat  (Fig. 
2),  a  substantial  part  of  the  thermal  heating  of  the  argon  pro¬ 
pellant  took  place  in  the  supersonic  part  of  the  nozzle,  which 
reduced  the  tnermal  thrust  contribution  (Fig.  6). 
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with  d^  and  dc  as  the  anode  and  cathode  diameter,  respec¬ 
tively.  The  maximum  thrust  TmfITiax  was  calculated  with  the 
largest  possible  anode  diameter  and  the  minimum  thrust  Tm>min 
with  the  smallest  possible  d^  (see  Fig.  2).  At  lower  currents, 
the  thrust  increases  almost  linearly  with  I  and  at  higher  cur¬ 
rents  attains  almost  the  slope  of  the  Tm  curve.  The  current 
and  thrust  are  limited  by  the  onset  condition. 

Onset  is  detected  visually  by  the  formation  of  bright  hopping 
spots  at  the  anode.  A  slight  increase  in  the  current  at  that 
point  leads  to  a  longer  attachment  of  the  spots,  causing  severe 
melting  wnich  destroys  the  anode.  In  the  experiments,  this  was 
avoided . 

The  specific  impulse  as  a  function  of  the  correlation  parameter 

l2/m  is  displayed  in  Fig.  12.  It  shows  almost  no  dependence  on 

the  mass  flow  and  reaches  a  maximum  of  about  11,000  m/s  at  the 

onset  points.  The  specific  impulse  is  in  a  wide  range  almost 

linear  with  I 2/m  and  parallel  to  a  line  between  I  and 

c  spm,max 

Ispm,min'  which  means  that  the  thermal  part  is  independent  of 
the  correlation  parameter. 


Fig.  12  Specific  impulse  as  a  function  of  the  correlation  pa¬ 
rameter  I  2/m  for  several  mass  flows  with  argon  ..as  propellant. 
ISpm,max  and  ^spn^min  are  the  calculated  magnetic  impulse 
for  maximal  and  minimal  anode  radius  (see  Fig.  2). 


In  Fig.  13,  the  arc  voltages  are  plotted  as  a  function  of  the 
discharge  current,  with  the  mass  flow  as  parameter.  Except  near 
onset,  the  arc  voltages  are  insensitive  to  mass  flow,  and  the 
increase  of  voltage  is  nearly  linear  with  the  current.  The 
start-up  voltage  is  about  25  to  30  V  at  1  kA.  The  onset  voltage 
of  roughly  70  V  is  nearly  the  sane  for  all  mass  flows  between 
0.3  and  0.8  g/s  and  in  the  range  of  the  measurements  in  Tokyo 
with  the  Kill  thruster13  at  0.6  g/s  argon. 
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Fig.  14  Thrust  efficiency  as  a  function  of  the^  specific  im- 
pulse  for  several  mass  flows  with  argon  as  propellant. 


With  the  measured  data,  a  thrust  efficiency  can  be  calculated  as 


2ml  *U 


This  thrust  efficiency  is  displayed  as  a  function  of  the  specif¬ 
ic  impulse  in  Fig.  14.  Because  the  Stuttgart  thruster  works  in 
the  predominantly  thermal  region,  the  thrust  and  therefore  the 
efficiency  rise  with  a  higher  propellant  flow. 
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E0  and  thermal  efficienc 


as  a  function  of  current 


Thermal  losses 

The  heat  loss  of  each  neutral  segment  and  of  the  electrodes  was 
measured  separately.  Fig.  15  shows  the  loss  factors,  defined  as 
heat  load  per  electrical  power  input  times  100,  of  anode  (V^) , 
cathode  ( Vc )  and  all  neutral  segments  (R^,  Rj,  throat  and  arc 
chamber),  eg  is  the  sum  of  all  thermal  loss  factors,  and  is 
connected  with  the  thermal  efficiency:  nxh  =  100  -  Eg , 
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wmch  is  also  plotted  in  this  figure.  At  low  currents,  the  loss 
factor  of  the  cathode  decreases  with  increasing  current  and  is 
almost  independent  of  current  near  the  onset  point.  The  loss 
factors  of  the  neutral  elements  show  only  a  small  increase  with 
rising  current,  even  near  the  onset  point,  which  similarly  indi¬ 
cates  only  a  small  temperature  rise  of  the  plasma  jet.  This 
means  that  the  severe  increasing  of  tne  anode  loss  factor  is 
almost  exclusively  due  to  an  increasing  anode  fall  voltage 
Ua*  Neglecting  radiation  and  possible  ablation  effects,  the 
energy  balance  of  the  anode  can  be  expressed  as:  iu 


VA  5  If 

=  U.  +  ^  -  T  +  $+  U 
I  A  2  e  e  conv , 


In  order  to  calculate  Ua  from  this  equation,  the  anode  heat 
loss  Qa  is  measured,  the  known  value  for  the  work  function  is 
used,  and  the  convective  heat  is  determined  from  the  measured 
heat  loss  of  the  neutral  segment  R2  next  to  the  anode.  From  in¬ 
vestigations  by  Maisenhalder  of  the  electron  temperature  in  the 
plasma  jet14  for  a  very  similar  thruster,  we  chose  an  electron 
temperature  near  the  anode  surface  of  15,000  K.  The  anode  fall 
voltage  Ua  is  shown  in  Fig.  16,  which  is  in  accordance  with 
tne  theoretical  thoughts  on  anodes  by  Hiigel.10  Tne  steep  rise  of 
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Fig.  16  C_a  leu  lated  anode _ fall  voltage  as  a  function  _of .  d  i  s- 

cnarge  current  for  0. 8  g/s  argon  as  propellant. 
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the  anode  fall  voltage  near  the  onset  point,  with  almost  un¬ 
changed  temperature  in  the  arc  jet,  points  to  a  changed  anode 
mechanism  at  onset,  which  could  also  be  detected  visually  during 
our  experiments. 

In  Fig.  17,  the  anode  loss  voltage,  defined  as  ^anode^1'  is 
shown  for  different  mass  flows  versus  the  correlation  parameter 
I2/m.  Again  the  steep  rise  of  the  anode  losses  near  the  onset 
point,  which  is  independent  of  the  mass  flow,  should  be  noted. 

The  thermal  efficiency  irrh  (Fig-  15)  of  this  thruster  is  about 
60%  and  decreases  near  the  onset  point  because  of  increasing  an¬ 
ode  fall  voltage,  as  discussed  before. 


Anode  gas  injection  (with  DT2  thruster) 


The  increase  in  the  anode  fall  voltage  near  onset  suggests  that 
this  phenomenon  is  at  least  partly  an  anode  problem.  Bakshtib 
and  Hiigel  1  u  built  their  onset  theories  thereon  as  an  anode  star¬ 
vation  process.  This  theory  and  experiments  at  Princeton16  sug¬ 
gest  that  gas  injection  near  the  anode  would  increase  the  criti¬ 
cal  I */m  to  higher  values.  The  anode  gas  is  injected  at  the  be¬ 
ginning  of  the  anode,  tangetial  to  the  anode  surface. 


Fig.  17  Anode  loss  voltage  as  a  function  of  the  correlation  pa 
raineter  I  Vifi  for  several  mass  flows  with  argon  as  propellant. 


voltage  of  3  V  was  taken  for  defining  onset  in  the  continuously 
running  MPD  thruster  (see  Fig.  18). 

In  order  to  justify  this,  the  current,  at  which  a  certain  anode 
loss  voltage  is  reached,  is  plotted  versus  the  injected  anode 
mass  fraction  with  the  anode  fall  voltage  as  parameter  (see 
Fig.  19). 

The  onset  rises  steadily  with  the  anode  gas  fraction.  The  arc 
voltage  (Fig.  20)  is  almost  the  same  at  lower  currents,  with 
somewhat  lower  values  for  the  higher  anode  gas  fractions.  At 
higher  currents,  the  voltage  curves  spread,  in  accordance  with 
the  higher  anode  fall  voltages.  This  increase  of  the  onset  cur¬ 
rent  must  concede  to  worse  thrust  values  at  higher  anode  gas 
fractions  (Fig.  21). 

To  overcome  this  handicap,  experiments  were  carried  out  with  the 
lighter  hydrogen  and  nitrogen  as  anode  gas,  where  the  masses 
were  chosen  so  that  the  particle  number  was  the  same  as  with 
argon  anode  gas.  Contrary  to  the  expectations  of  the  anode  star¬ 
vation  theories,  these  anode  gases  changed  neither  the  voltage 
nor  the  onset  behavior,  and  showed  results  similar  to  those  in 


experiments  without  any  anode  gas  injection.  This  phenomenon  can 
only  then  be  explained  by  Hiigel's  theory  if  the  injected  anode 
gases  are  not  ionized  at  all.  This  requires  a  very  low  ioniza¬ 
tion  rate  coefficient  of  the  injected  gases. 

In  Figs.  22  through  25  the  voltage  current  characteristics  are 
drawn  in  a  double  logarithmic  manner  for  argon  as  main  (cathode) 
propellant  and  with  different  elements  as  anode  gas.  The  total 
mass  flow  in  each  case  is  0.8  g/s.  The  curves  all  show  a  similar 
oenav ior--they  rise  linearly  to  a  certain  point  of  about  4.1  kA, 
then  change  abruptly  to  a  steeper,  also  constant,  slope. 

The  first  part  up  to  the  angle  point,  which  indicates  fully  ion¬ 
ized  plasma,  is  practically  identical  and  is  independent  of  the 
anode  gas  fraction.  The  second,  steeper  part  of  the  voltage 
characteristics  depend  on  the  anode  gas — the  more  anode  gas,  the 
less  steep  is  the  slope.  An  explanation  may  be  that  the  injected 
cold  gas  must  be  ionized  and  so  the  voltage  rise  is  smaller  than 
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Fig.  21  Thrust  vs.  discharge  current  with  different  portions  of 
anode  gas,  tor  a  total  mass  flow  of  0.8  g/s  argon . 
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Fig.  22  Arc  voltage  vs.  discharge  current  with  different  por- 
s  of  argon  anode  gas,  for  a  total  mass  flow  of  0.8  q/s 
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Fig.  23  Arc  voltage  vs.  discharge  current  with  different 


tions  of  nitrogen  anode 


for  a  total  mass  flow  of  0.8 
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Fig.  24  Arc  voltage  vs.  discharge  current  with  different  por¬ 
tions  of  hydrogen  anode  gas,  for  a  total  mass  flow  of  0.8  g/s. 
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Fig.  25  Arc  voltage  vs.  discharge  current  with  xenon  anode  gas 
for  a  total  mass  flow  of  0.8  q/s. 
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without  ionization.  In  contrast  to  this  fact,  with  smaller  atom¬ 
ic  masses  (that  means  higher  particle  numbers  must  be  ionized  at 
constant  mass  flow)  the  slope  becomes  steeper.  This  eventually 
can  be  attributed  to  the  better  heat  conduction  with  gases  of 
low  atomic  mass. 

Without  anode  gas,  the  onset  was  observed  visually,  as  described 
above,  by  detecting  arc  attachment  spots  at  the  anode.  With  an¬ 
ode  gas  injection,  the  rise  of  the  anode  fall  voltage  is  driven 
to  higher  currents  and  the  onset  is  split:  before  the  appear¬ 

ance  of  the  anode  spots,  the  discharge  begins  to  oscillate  with 
a  frequency  of  about  220  kHz.  With  a  further  rise  in  the  cur¬ 
rent,  the  frequency  rises  and  becomes  more  stochastic.  For  exam¬ 
ple,  witn  a  5%  anode  gas  fraction: 

Begin  of  oscillations:  I  =  4350  A,  f  =  222  kHz 

at  I  =  4540  A,  f  =  267  kHz 

and  I  =  4700  A,  f  =  310  kHz. 

In  Fig.  26  an  oscillogram  of  such  an  instability  is  shown,  taken 

with  a  20  MHz  sample  rate  transient  recorder  from  signals  of  a 
floating  potential  probe  in  the  plasma.  This  split  of  the  onset 
(with  the  pulsed  thrusters  there  even  exists  another  onset 
point,  namely  the  angle  point  of  the  voltage  current  character¬ 
istic  as  shown  in  section  III)  indicates  that  the  onset  is  not 
only  an  anode  phenomenon,  but  also  a  stability  problem,  as  ex¬ 
plained  by  a  theory  of  Schrade  at  the  IRS.X/  The  beginning  of 
the  oscillation  is  a  function  of  the  anode  gas  fraction  as  shown 
for  the  case  Ar/Ar  (cathode  gas/anode  gas),  mtot  =  0*8  g/s: 


100/  0 
98/  2 
95/  5 
90/10 

Unfortunately,  no  performance  measurements  in  the  unstable  re¬ 
gion  could  be  obtained,  because  the  thrust  and  thermocouple 
measurements  are  disturbed  by  these  oscillations,  and  a  remedy 
has  not  yet  been  found. 
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mass  flow.  With  higher  discharge  currents,  the  thrust  values  ap¬ 
proach  each  other.  The  efficiencies,  shown  in  Fig.  29,  are  plot¬ 
ted  against  the  specific  impulse.  The  lower  efficiencies  of  N  i 
are  due  to  the  higher  voltage,  and  the  decline  with  higher  I3p 
corresponds  to  the  steep  rise  of  tne  voltage  near  onset. 

With  anode  gas,  the  results  are  in  general  similar,  as  demon¬ 
strated  in  Figs.  30  to  32.  Again  the  shifting  of  the  onset  to 
higher  current  values  can  be  observed,  where  it  is  interesting 
to  note  that  the  slope  after  the  angle  point  with  N  ^/N ;  is  re¬ 
markably  smaller  than  with  Ar/Ar  or  with  the  combination  Ar/N  z 
in  Fig.  23. 
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Fig.  .19  Thrust  efficiency  as  a  function  of  the.,  specific  ..im¬ 
pulse  lor  an  egual  volume  flow  of  argon  (0.8  g/s)  and  nitrogen 
(0.S6  j,  sj  as  propellant,  without  anode  gas. 
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Fig.  31  Measured  thrust  vs.  current  for  an  equal  vol ume  flow 
of  a  r  jon__  ( U_.  8  g/s]  and  nitrogen  (0.56  g/s)  with  an  anode  gas 
mass  tract  ion  of  1iU.  Tm>(nax  and  Tmfmin  are  the  calculated 
electromagnetic  thrust  for  the  maximal  and  minimal  possible 

anode  radii. 
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Fig.  32  Thrust  efficiency  as  a  function  of  the  specific  impulse 
for  an  equal  volume  flow  of  argon  (0.8  g/s)  and  nitrogen  (0.56 
g/  sj  as  propellant,  with  an  anode  gas  mass  fraction  of  10%. 


2 . 2  Cylindrical  MPD  Thruster  with  Elongated  Electrodes 

In  contrast  to  the  nozzle-type  MPD  thrusters,  which  are  in  ef¬ 
fect  Hybrids  of  pure  MPD  and  thermal  arcjets,  cylindrical 
thrusters  are  predominantly  MPD  devices.  Further,  the  elongated 
electrodes  are  hoped  to  reduce  the  pinch  effect  found  to  be  an 
essential  factor  affecting  the  anode  mass  starvation,  according 
to  tne  calculations  of  Hugelil"  and  the  Princeton  experiments.1 
Tnese  catnodes  wiLl  be  mostly  radiation-cooled  and  may  be  more 
limited  in  our r en t - rat  r  y  l  ng  capability.  The  tntuster  has  been 
ies  igned  based  on  available  information  from  previous  experi¬ 
ments  '  i  n  o  i  j.lmj  Princeton  data)  and  calculations  tor  maximum 

I'  m  and  .specif  ic  impulse,  in  tins  first  version  of  a  continu¬ 
ously  running  cylindrical  thruster  with  a  long  cathode,  the  an¬ 
ode  an i  all  neutral  segments  will  be  water-cooled  ' Fig.  13).  A 
;»*cond  ,  non-cooled  version  wi»r.  t  fie  same  electrodes  will  be 
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It  is  not  yet  evident  whether  the  full  length  cathode  will  oper¬ 
ate  smoothly  in  a  steady-state  mode.  Therefore,  the  cathode  is 
being  built  as  a  displaceable  one,  where  the  discharge  will  be 
ignited  in  a  retracted  position.  When  the  cathode  tip  begins  to 
glow,  the  cathode  will  be  protruded  to  the  foreseen  position. 

The  anode  is  divided  into  three  parts  to  give  some  insight  into 
the  current  density  distribution  on  the  anode  and  to  compare  it 
with  the  theoretically  evaluated  values. 

Since  the  anode  rings  can  be  switched  separately,  this  thruster 
with  displaceable  cathode  will  allow  a  wide  range  of  geometric 
variations  to  ne  tested.  Experiments  at  Princeton^u  showed 
strong  dependencies  of  tne  onset  1 2 /m  values  on  the  place  where 
the  propellant  was  fed.  With  this  thruster  design,  therefore, 
two  propellant  ports  are  anticipated  in  a  slit  around  the  cath¬ 
ode  and  in  a  slit  at  the  anode  radius,  both  at  the  backplate.  A 
variation  of  the  propellant  fractions  at  these  ports  should  show 
whether  the  observed  onset  variations  with  quasi-steady  pulsed 
configurations  can  be  reproduced  with  a  steady-state  device. 

In  Figs.  34  and  35  first  results  of  the  influence  of  the  cath¬ 
ode  position  variation  on  the  current  distribution  are  shown. 
The  propellant  flow  was  a  total  of  2  g/s  argon,  divided  3:1  re¬ 
spectively  between  the  cathode  and  the  outer  (anode)  injection 


Fig.  it  Cylindrical  MFD  thruster  head  schematic 


port.  In  Doth  cases,  the  arc  was  ignited  at  a  position  of  ca .  -5 
mm  (0  =  backplate)  with  about  1200  A;  the  cathode  was  shifted  to 
ca .  +5  mm  and  there  the  current  was  increased  to  the  planned 
current  level,  2000  and  3000  A.  Then  the  cathode  was  fed  by 
motor  with  a  velocity  of  *  1  mm/s  downwards.  Both  cases  show 
similar  behav ior-- f i rs t  the  rearmost  anode  ring  took  almost  all 
the  current,  the  middle  and  front  ring  only  very  little.  When 
the  cathode  tip  reached  the  region  of  the  rear  anode  ring,  the 
anode  current  distribution  changed  dramat  ically--the  main  cur¬ 
rent  part  switched  from  the  back  ring  to  the  front  ring,  while 
the  middle  anode  ring  was  almost  passed  over.  Both  tests  had  to 
oe  stopped  when  the  cathode  tip  reached  the  beginning  of  the 
first  anode  ring  ca .  30  mm  behind  the  front  plate,  and  relative¬ 
ly  weak  anode  attachment  spots  occured  at  the  outer  edge  of 
anode  ring  1.  The  voltage  first  decreased  slightly  witn  the  pro¬ 
truding  of  the  cathode  and  remained  then  practically  constant. 
With  the  lower  current  (2000  A)  at  15  mm  a  rearrangement  of  the 
discharge  was  ooserved,  accompanied  by  a  voltage  increase.  Dur¬ 
ing  ignition,  some  backward  flashing  of  the  arc  into  the  chamber 
in  front  of  the  throat  occured,  which  caused  some  damage.  As 
these  first  experiments  should  only  verify  the  working  of  a  con¬ 
tinuously  running  cylindrical  MPD  arc,  no  thrust  was  measured. 
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Fig.  37  Nozzle-type  MPD  thruster  with  rad xat ion-cooled  tungsten 
anode  (HAT ) . 

2.3  Hot  Anode  Thruster 


A  substantial  part  ( *  30%)  of  the  input  power  of  an  MPD  thrust¬ 
er  is  consumed  in  the  anode.  Thus,  a  decrease  in  these  anode 
losses  would  heighten  the  performance  of  these  devices. 

MPD  tnrusters  with  radiat ion-cooled  anodes  show  higher  thermal 
°t  f iriencies  than  those  with  water-cooled  ones,  but  they  were 
used  only  with  thrusters  of  a  much  lower  power  level  and  mainly 
of  the  applied  magnetic  field  type.  b,  '  ‘ 

Tne  "Hot  Anode  Thruster”  (HAT)  is  an  attempt  to  use  these  advan¬ 
tages  of  rad  lat  ion-cooled  anodes  tor  thrusters  in  the  power 
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class  of  200  kW.  Its  schematic  is  shown  in  Fig.  17.  In  order  to 
provide  a  sufficiently  large  area  to  radiate  the  anode  input 
power,  the  anode  flange  has  a  diameter  of  roughly  40  cm.  It  is 
made  out  of  tungsten;  the  areas  not  intended  to  be  loaded  by  the 
current  are  covered  by  tantalum  carbide,  which  is  a  relatively 
weak  insulator.  This  coating  insulates  the  radiating  anode  areas 
and  at  the  same  time  increases  the  coefficient  of  emmissivity  by 
about  a  factor  of  2  to  3  versus  pure  tungsten. 

In  a  first  prototype,  the  hot  anode  is  shielded  from  the  other 
parts  of  the  thruster  by  a  water-cooled  copper  plate,  which  is 
diffusion  welded  and  forms  one  of  the  neutral  rings  of  the 
nozzle . 

In  order  to  compare  these  experimental  results  with  those  of  a 
water-cooled  one,  the  thruster  is  designed  as  closely  as  possi¬ 
ble  in  the  form  of  the  "normal"  Stuttgart  nozzle-type  thruster. 
The  current  is  fed  to  the  anode  by  16  screws  made  of  TZM ,  a 
molybdenum  alloy,  which  at  the  same  time  act  as  a  mounting  sup¬ 
port  for  the  heavy  anode  flange. 

In  Fig.  38,  the  result  of  a  thermal  analysis  is  shown  for  the 
case  of  60  kW  power  input,  which  is  assumed  to  load  the  anode  in 
a  linear  manner.  The  analysis  was  performed  using  a  finite  ele¬ 
ment  proced ure , 2  1  '  2  J  which  utilizes  the  SMART  temperature  calcu¬ 
lation  code,  developed  at  the  Institut  fur  Statik  und  Dynamik 
(ISD)  of  the  University  of  Stuttgart.  As  a  result,  it  is  shown 
here  that  even  at  such  a  high  power  input  (which  corresponds  to 
roughly  200  kW  electrical  power  input  to  the  thruster)  the  tem¬ 
perature  would  not  at  any  place  exceed  2600°  C. 

The  initial  tests  with  the  HAT  are  intended  to  be  carried  out 
in  a  newly  installed  stainless  steel  vacuum  tank  of  2  m  in  diam¬ 
eter  and  6  m  in  length.  It  is  double  walled  to  provide  cooling 
capab  i  I l ty . 


3 .  Quasi-Steady  Operated  MPD  Thrusters 


In  Fig.  39,  specific  impulse  data  from  our  experiments  (station¬ 
ary  mode)  and  from  some  Princeton",  and  Tokyo  1  6  quasi-steady 
experiments  were  plotted  against  the  correlation  coefficient 
I  /m .  The  main  differences  that  are  striking  are  the  substanti¬ 
ally  smaller  I Vm  values  and  therefore  the  smaller  specific  im¬ 
pulses  which  are  reached  in  the  continuous  mode. 

To  check  whether  these  great  differences  between  these  pulsed 
thrusters  and  our  continuous  thrusters  are  a  function  of  the  de¬ 
sign  configuration  or  of  the  operation  mode,  thrusters  identical 
in  geometry  and  electrode  material  were  built  (Figs.  40  and  41) 
and  investigated  in  the  quasi-steady  mode . 2  3  The  neutral  seg¬ 
ments  which  form  the  nozzle  (Figs.  2  and  3)  were  replaced  by 
ceramic  insulators  (Macor).  The  power  source  was  a  capacitor 
bank  with  electrolytic  capaci tors  1  b  of  approximately  6  kJ,  which 
provides  a  quas  i-stationary  discharge  period  for  about  2  ms. 
Since  the  vacuum  system  was  sufficiently  large,  the  gas  flow  did 


Fig.  39  Specific  impulse  as  a  function  of  the  correlation  pa- 
rameter  I  1 /m  for  different  thrusters,  stationary  and_  pulsed 

ones : Stuttgart  thruster,  stationary  (measured  thrust)  and 

pulsed  (calculated  thrust),  Tokyo  K  III  thruster,  Princeton 
benchmark  and  extended  anode  thruster.  Propellant  is  in  all 
cases  argon. 


not  nave  to  be  switched;  tnat  means  that  the  electrical  pulse 
was  "switched  into"  a  stationary  yas  flow.  The  main  discharge 
was  triggered  by  a  high  voltage  spark,  which  was  ignited  between 
the  cathode  and  the  copper  arc  chamber  wall.  The  discharge  cur¬ 
rent  was  measured  by  a  Pearson  current  transformer  (301X)  and 
the  discharge  voltage  at  the  thruster  terminals  with  a  fast  iso¬ 
lation  amplifier.  Botn  signals  were  stored  on  transient  record¬ 
ers.  The  main  and  the  ignition  power  supply  potentials  were 
floating  . 


Cathode 


Arc  Chamber  (Copper) 


Fig.  4  0  Thruster  head  configuration  for  quasi- 
cor respond  ing  to  thruster  DTK 
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Fig.  4  1  Thruster  head  configuration  f  or__cjuas  i-s  teady  operation , 
cor  respond  mu  to  steady-state  thruster  DT3. 


As  onset  criterion,  the  appearance  of  nigh  voltage  oscillations 
on  the  voltaye  traces,  which  happen  to  coincide  with  the  usual 
10%  fluctuation,  was  taken.  In  Fig.  42,  the  result  of  the  meas¬ 
ured  relative  voltaye  fluctuations  for  different  mass  flows  are 
plotted  versus  the  correlation  parameter  I‘/m.  For  1  /m  <  1*10 
A"s/kg,  no  voltage  oscillations  could  be  observed,  but  they 
start  almost  instantaneously  above  this  "critical"  value  with  a 
magnitude  of  10%  and  above,  so  that  for  this  thruster  and  an  am¬ 
bient  pressure  of  ca .  2  *  1 0 —  mbar,  the  critical  or  onset  point 

proved  to  be  1  •  1 0 1 0  A"s/kg,  a  factor  of  2.5  smaller  than  within 
the  continuous  operation  mode! 


Fig.  42  Relative  voltage  fluctuations  for  different  mass  flows 
v s.  correlation  parameter  I2/m,  with  argon  as  propel lant . 


In  Figs.  43  to  45,  typical  current  and  voltage  wave  forms  are 
plotted  from  the  transient  recorder.  Fig.  43  shows  the  traces 
for  a  typical  discharge  under  the  critical  region  ( I  Vm ”0 . 5 • 1 0 1 u 
A"s/kg):  the  current  pulse  is  almost  rectangular  and  the  volt¬ 
age  pulse  is  also  smooth  and  over  a  long  period  constant. 

Fig.  44  shows  the  current  pulse  and  two  voltage  pulses  for  the 
region  just  above  the  onset  (1.0*10iU  <  I  2/m  <  1.5*10lLl  A^s/kg): 
the  current  pulse  deforms  and  the  constant  part  gets  smaller.  In 
Fig.  45,  the  case  high  above  onset  is  plotted:  the  current  shows 
a  distinct  peak  at  the  beginning  which  corresponds  to  the  de¬ 
pression  in  the  voltage  trace.  The  whole  pulse  period  is  pro¬ 
longed  . 

Current  voltage  characteristics  with  the  mass  flow  as  parameter 
are  shown  in  Fig.  26  in  a  double-log  plot.  They  show  the  famil¬ 
iar  kinked  form13,2'  with  a  constant  slope  up  to  a  certain  cur¬ 
rent  value  and  then  an  abrupt  jump  to  a  steeper,  also  constant 


b)  c) 


tig.  44  Typical  a)  current  and  b) ,  c)  voltage  wave  forms  just 
above  onset  .  b)  I  2/m  =  1.1  *10  10  A2s/kg,  c)  I  2/m  =  1.4  *10  lu 
A  1 s/ kg  )  . 


Fig.  4  5  Typical  a)  current  and  b)  voltage  wave  forms  well 
above  onset  (  I  Vm  =  2.9  *10  10  A‘s/kg). 


Current  during  constant  Phase  [  kA  ] 
for  different  Mass  Flows 

Fig.  46  Qirrrent-  voltage _ character's  t  ics  with  the  argon  mass 

flow  as  parameter.  Dotted  line  represents  the  continuous  thrust¬ 
er  as  comparison. 


slope.  In  tne.se  experiments,  this  knee  coincides  with  the  on¬ 
set  current,  wriicli  seems  reasonable,  since  the  “onset"  marks 
another  physical  n-luviot  ot  tne  arc  with  higher  anode  losses 
and  voltages.  Kuriki'  called  this  knee  "critical  carrent",  bat 
defined  as  onset  a  "limitin'.)  current"  which  is  located  somewhere 
on  toe  steep  branch  ot  tne  character  i  st  ic  .  Choueiri  1  found  in 
tne  "knee"  tne  beginning  ot  an  instability  which  is  related  to 
the  transition  to  a  fully  ionized  plasma,  and  which  is  damped 
away  at  rusher  currents  still  well  below  onset  point.  As  in  the 
case  with  steady  state  thrusters,  our  experiments  with  these 
pulsed  devices  showed  that  currents  above  onset  point  are  always 
connected  with  melted  spots  at  the  anode,  because  the  former 
diffuse  anode  attachment  of  the  arc  probably  constricts  to  over¬ 
come  a  so-called  "anode  starvation".  10  We  have  noticed  such 
melt  in*)  spots  wi*-.n  all  currents  above  tne  "knee”,  so  tnat  set¬ 
ting  tnis  "knee  current"  as  onset  current  seems  meaningful. 

In  the  same  figure  (4b),  a  dotted  line  is  used  to  plot  the  volt¬ 
age  characteristic  of  the  continuous  mode;  this  dotted  line  has 
lower  values  and  the  oscillations  appear  after  this  angle 
point..  This  mdicates  that  these  two  instabilities  may  be  ot 
different  origin.  Tne  lower  voltages  are  probably  caused  by  the 
lower  cathode  fall  with  a  glowing  tungsten  cathode  compared  to  a 
cold  one.  A  furtner  indication  that  the  onset  point  is  a  change 
in  the  physical  behavior  is  shown  in  Fig.  47.  The  current  output 
of  tne  PFN  for  a  given  charging  voltage  U  }  as  a  f met  ion  of  the 
mass  flow  shows  a  dist  met  kmk  at  a  current,  which  coincides 
with  the  onset. 


In  a  series  of  experiments,  the  dependence  of  onset  with  the  am- 
b  mu!  pt  es'iut  e  was  invest  igated ,  which  showed  that  with  increas¬ 
ing  nmnient  pressure  the  onset  current  rises  substantially 
(Fig.  4b  i  .  This  result  contradicts  the  findings  with  steady- 
state  thrusters,  wmen  indicates  that  both  operation  modes, 
qu as i- steady  pulsed  and  continuous,  are  not  only  distinguished 
by  differe-i*  -at  hod*'  conditions,  i.e.  cold  and  hot  ,  but  also  by 
o  t  n >•  r  Oil!  m  Known  •■fleets. 
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Thrust 

Thrust  is  measured  in  an  indirect  manner.  The  thruster  is  mount¬ 
ed  on  a  pendulum  tnrust  oalancu  and  the  total  impulse  of  the 
discharge  is  measured  by  the  pendulum  displacement.  The  thrust 
is  then  calculated  by  assuming  a  rectangular  discharge  pulse. /0 
The  balance  is  calibrated  by  applying  a  known  impulse. 

Measurements  of  thrust  were  taken  with  argon  (Fig.  49)  and  ni¬ 
trogen  (Fig.  50)  as  propellant  for  several  mass  flow  rates  and 
currents.  In  both  figures  the  thrust  of  the  continuously  running 
thruster  is  displayed  for  comparison  as  dashed  lines  for  one 
mass  flow,  along  witn  the  possible  minimal  and  maximal  magnetic 
tnrust.  The  thrust  of  the  pulsed  experiments  is  smaller,  espe¬ 
cially  at  lower  currents.  A  possible  explanation  could  be  that 
the  thermal  heating  of  the  propellant  is  incomplete  because  of 
turbulences  caused  by  pumping  cathode  attachment  spots  of  the 
arc.  Tms  is  confirmed  by  high  speed  films  (20000  f rames/s) --no 
stabilization  of  the  jet  plume  could  be  observed  during  the  2  ms 
pulse.  The  specific  impulse  as  a  function  of  the  parameter  l2/m 
is  plotted  in  Fig.  51  with  mass  flows  as  parameter  for  nitro¬ 
gen.  The  specific  impulse  runs  almost  linear  to  I vm  over  a  wide 
range  and  runs  parallel  to  the  magnetic  specific  impulse. 

The  efficiencies,  shown  in  Figs.  52  and  53,  are  lower  than  in 
the  continuous  case  because  of  the  lower  thrust  and  the  higher 
discharge  voltages. 


Thrust  as  a  function  of  current  with  nitrogen  mass  flow 
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Magnetic  Probe  Measurements 


Magnetic  probe  measurements  are  not  possible  on  a  continuous 
plasma  jet  because  of  the  high  heat  loads  on  the  probe.  These, 
however,  are  of  minor  influence  with  pulsed  jets.  With  an  iso¬ 
lated  magnetic  probe,  the  current  contour  lines  were  measured 
and  are  compared  with  calculated  lines  in  Fig.  54.  They  coincide 
fairly  well  for  the  10%  and  20%  current  contour  line.  The  fact 
that  in  the  measured  case  the  anode  attachment  is  more  concen¬ 


trated  downstream  of  the  nozzle  end  than  in  the  calculated  case 


may  be  attributed  to  the  theoretical  model  which  neglects  the 
boundary  layer  along  the  nozzle  and  anode  walls  and  taken  a  uni¬ 
form  conductivity  over  the  cross  section.  The  theoretical  model 
is  described  in  detail  in  chapter  5. 
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Fig.  54  a)  Calculated  and  b)  measured  current  contour  j^ines  of 
an  MPD  discharge  in  the  nozzle-type  plasma  thruster  D2  for  0.8 
g/s  argon  and  1  300  A  discharge  current . 
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With  the  measured  current  contour  lines  and  a  calculated  elec¬ 
trical  conductivity,  taken  as  constant,  the  Joule's  heat  source 
distribution  was  determined  in  Fig.  30  a)  for  1300  A  and  in  b) 
for  3150  A.  With  higher  current,  the  heating  zone  constricts  to 
the  middle  axes. 
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Fig.  5  5  Isol  i  nes  for  the  Ohmic  heat  source  term  for  a)  1300  A 
and  b)  3150  A  discharqe  current. 
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Concluding  Remarks 

The  experimental  results  presented  are  part  of  a  greater  effort 
to  understand  MPD  propulsion  and  the  related  phenomena,  which 
include  theoretical  investigations,  both  of  the  midstream  flow 
field  and  electrodes  and  of  erosion.  The  goal  of  the  experiments 
is  twofold:  first  to  give  a  basis  of  data  for  the  theories  and 
second  to  come  up  with  better  and  more  efficient  thruster  de¬ 
signs,  since  the  performances  of  self-field  MPD  thrusters,  at 
least  in  the  investigated  power  field  <  500  kW,  are  quite  mod¬ 
est.  With  regard  to  this  second  goal,  the  experiments  presented 
here  can  be  summarized: 

1.  The  acceleration  machanism  of  nozzle-type  MPD  thrusters  in 
the  tested  performance  regime  is  not  only  determined  by  the 
self  magnetic  forces  but  also  to  a  great  extent  by  thermal 
expansion . 

2.  The  influence  of  cathode  position  on  the  thruster  performance 
is  a  further  proof  of  the  important  role  of  the  thermal  ac¬ 
celeration  . 

3.  A  variation  of  the  ambient  tank  pressure  up  to  0.5  mbar  hard¬ 
ly  influences  the  thruster  performance.  Also,  the  onset  point 
is  not  shifted  with  steady  operation  mode.  Thus,  the  evalua¬ 
tion  of  specific  impulse  and  efficiency  with  direct  thrust 
measurement  are  valid;  no  entrainment  of  ambient  gas  is  to  be 
expected  at  operation  below  onset. 

4.  The  performance  improvement  hoped  for  in  nozzle-type  thrust¬ 
ers  with  anode  gas  injection  due  to  shifting  of  the  onset  to 
higher  currents  was  not  realized. 

5.  With  anode  gas  injection,  the  onset  phenomena  were  split: 
at  least  two  onset  effects  with  steady  state  MPD  thrusters 
could  be  observed: 

a)  anode  spot  attachment  due  to  particle  starvation 

b)  high  frequency  current  oscillations. 

This  indicates  that  the  onset  limitations  are  caused  by  sev¬ 
eral  different  physical  reasons. 

6.  At  these  low  power  ranges,  the  lighter  propellant  nitrogen 
snows  no  advantages  over  argon. 
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7.  First  experiments  with  a  cylindrical  thruster  proved  the 
feasibility  of  a  cathode  protruding  from  the  backplate.  The 
current  distribution  on  the  anode  depends  to  a  large  extent 
on  the  cathode  position. 

8.  Quasi-steady,  pulsed,  nozzle-type  thrusters  showed  worse  per¬ 
formances  in  the  power  range  considered  than  steady-state  op¬ 
erated  thrusters: 

a)  higher  voltages,  due  to  the  cold  cathode 

b)  lower  thrust,  perhaps  due  to  the  worse  input  of  the  ohmic 
heat  because  of  turbulent  flow  in  the  thruster 

c)  lower  onset  values,  which  occur  at  the  "knee"  of  the  volt¬ 
age  current  characteristics. 
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4.  Cathode  Erosion  Studies  on  Ml>u  Thrusters 

4  .  1  Introduct ion 

It  is  a  known  phenomenon  that  in  a  pulsed  MPD  thruster  1 ^  iJ,  11 
and  during  the  starting  phase  in  a  continuous  MPD  thr us  ter  i  2 ,  i i 
the  arc  cathode  attachment  resembles  many  hot  spots  of  high  cur¬ 
rent  density.  These  spots  are  usually  highly  nonstationary  and 
jump  across  the  cathode  surface.  Such  a  spotty  attachment  is  in¬ 
deed  the  most  common  one  on  cold  metal  cathodes  of  low  pressure 
arcs'34,  35  and  seems  to  be  an  intrinsic  feature  of  those  dis¬ 
charges.  Despite  the  growing  literature  on  this  phenomenon,  it 
is  still  not  fuLly  understood.  However,  since  the  time  resolu¬ 
tion  of  electrical  and  optical  measuring  techniques  went  up  into 
the  ns  and  even  beyond  the  nanosecond  range,  and  since  scanning 
electron  microscopy  became  a  widespread  technique  of  investi¬ 
gating  the  arc  traces  left  on  cathode  surfaces,  one  has  learned 
much  more  on  the  mechanisms  of  this  spotty  attachment.  Today  one 
distinguishes  between  a)  microspots  with  diameters  of  several  iim 
(maybe  even  less)  and  with  current  densities  of  the  order  of 
10 12  A/m 2  (maybe  even  higher)  and  displacement  velocities  of 
several  hundred  meters  per  second  and  b)  spots  which  are  much 
more  stationary  (velocities  of  the  order  of  several  m/s)  with 
diameters  of  up  to  100  iim.Jb  These  latter,  larger  spots  appear 
to  be  the  result  of  spot  clustering;  i.e.,  due  to  favorable 
emission  conditions,  the  microspots  do  not  spread  over  a  larger 
area  but  rather  stick  to  a  smaller  area,  and  by  locally  over¬ 
heating  this  area  cause  larger  molten  zones  and  relatively  heavy 
damage  to  the  cathode.  If  such  a  clustering  does  not  occur,  the 
inicrospot  damage  is  less  severe  because  of  the  rapid  spot  mo¬ 
tion,  and  because  consequently  there  is  more  widespread  averaged 
heatload  on  the  cathode  surface. 

Another  quite  interesting  feature  is  the  fact  that  the  average 
number  of  spots  is  proportional  to  the  arc  current,  and  that  a 
single  spot  can  carry  only  a  limited  current,  the  quantity  of 
which  depends  on  the  cathode  material.  1  ' ,  !  4 


Within  a  transverse  magnetic  tield,  these  spots  jump  preferably 
into  the  retrograde  direction,  i.e.  they  move  oppositely  to  Am- 

j  9  4  o 

pere's  law  as  long  as  the  ambient  pressure  is  small  enough.  , 

These  are  only  some  of  the  most  striking  known  and  established 
effects  observed  on  spotty  cathodes.  Many  more  detailed  spot 
features  have  been  observed  and  put  forward,  but  they  cannot  be 
dealt  with  in  the  framework  of  this  paper.  Reference  shall 
therefore  be  made  to  the  literature.35 


Nevertheless,  since  a  variety  of  these  effects  are  obviously  in¬ 
volved  at  the  cathode-arc  interphase  of  MPD  thrusters,  it  is 
mandatory  that  this  spot  behavior  be  investigated  in  order  to 
understand  and  predict  the  discharge  and  erosion  effects  during 
the  starting  phase  of  a  continuous  or  during  the  pulse  phase  of 
a  quasi-steady  MPD-arc  thruster. 

After  the  starting  phase  of  a  continuously  running  thruster, 
which  lasts  about  a  fraction  of  a  second  or  a  few  seconds  de¬ 
pending  on  the  arc  current,  the  cathode,  which  usually  consists 
of  a  thoriated  tungsten  rod — so  far  the  only  successfully  used 
in  continuous  thrusters--starts  to  glow,  and  the  original  spotty 
attachment  becomes  a  quiet,  more  diffuse  one.  The  current  densi¬ 
ties  of  these  diffuse  attachments  on  bright  glowing  tungsten 
cathodes  are  of  the  order  of  10  /  A/m 2  and  are  commonly  explained 
by  thermionic  emission,32  while  the  high  current  densities  of 
the  spots  can  only  be  interpreted  by  means  of  field  or  enhanced 
thermal  field  emission.41,42 


The  following  section  reports  about  erosion  measurements  on 
pul  seii  ano  rout  muons  MPD  thrusters  and  consequently  on  spotty 
and  diffuse  cathode  attachments,  while  the  subsequent  sections 
give  a  partially  quantitative  theoretical  interpretation  and 
discussion  of  the  experimental  findings. 


The  cathode  erosion  tests  have  been  conducted  on  the  nozzle-type 
MPD  thruster  DTI,  shown  in  Fig.  2.  For  the  pulsed  thruster  ex¬ 
periments,  the  water-cooled  segments  up  to  the  anode  ring  have 
been  replaced  by  a  heat  resistant  glass  ceramic  (MACOR)  and  the 
water  cooling  in  anode  and  cathode  not  applied.  The  eroded  mass 
is  determined  by  weighing  the  cathode  before  and  after  the  ex¬ 
periment  with  a  precision  balance  having  an  accuracy  of  1  mg  at 
a  weighing  range  of  1  kg. 

Steady  State  MPD  Thruster 


With  steady  state  operating  MPD  tnrusters,  two  different  modes 
of  erosion  must  be  distinguished: 

a)  Starting  with  a  cold  cathode,  heavy  erosion  takes  place, 
which  can  usually  be  detected  as  a  sort  of  "spitting"  of  glowing 
particles  or  droplets.  Some  tenth  of  a  second  to  several  seconds 
later,  the  cathode  tip  becomes  hot  and  glowing  so  that  thermion¬ 
ic  or  enhanced  thermal  field  emission  can  occur  and  the  dis¬ 
charge  becomes  stable  and  quiet. 

b)  During  the  steady  state  operation  phase,  the  arc  attachment 
is  confined  to  the  bright  glowing  area  of  the  cathode  tip.  The 
erosion  effect  is  much  smaller,  and  no  glowing  particles  are 
ejected . 

To  measure  these  different  erosion  rates,  two  series  of  tests 
were  performed:  long  duration  and  start-only  runs. 


Erosion  during  starting  phase.  During  the  "start-up  phase",  the 
arc  attachment  consists  of  many  highly  unstationary  hot  spots 
which  jump  over  the  cathode  surface  and  which  cause  a  relatively 
high  evaporation  and  even  splashing  of  cathode  material. 

A  current  voltage  plot  of  this  initial  starting  phase  (0  t  < 
t  x  )  is  shown  in  Fig.  56.  To  determine  the  average  eroded  mass 
during  this  starting  phase,  experiments  with  30  to  50  runs  with 
different  preset  ignition  currents  were  performed.  In  each  case, 
the  arc  was  extinguished  after  stabilization.  To  exclude  eventu¬ 
al  thermal  effects  in  the  starting  phase,  a  waiting  period  of 


I  n  K 1 '  i  • .  •>•*{  and  >'* ,  •  uvet  a  g*-  mas  loss,  is  plotted  alter  about 
10  to  ad  r  1:1 ,  as  a  t  a  not  ion  of  the  ign  it  ton  current  and  the  cor¬ 
respond  inj  acc  u  ii  a  1  u 1  ed  charge  Jar  inj  the  starting  phases,  re¬ 
spect  ively .  Trie  total  average  mass  loss  per  arc  ignition  amounts 
to  about  2a  mg  for  a  preset  arc  ignition  current  of  600  A  and 

decreases  to  about  10  mg  and  to  a  little  less  if  the  ignition 

current  is  1000  A  and  more.  Since  with  increasing  current  lj 
tne  cathode  heats  up  faster,  and  therefore  the  time  elapsed  dur¬ 
ing  the  starting  phase  becomes  smaller,  the  total  mass  loss  de¬ 
creases  with  increasing  ignition  current.  As  a  function  of  the 
accumulated  electric  charge  Qj  ,  however,  the  mass  loss  in¬ 
creases  witn  Of  almost  Linearly.  By  taking  the  least  square 
fit  slope  wo  obtain  an  average  erosion  rate  for  the  starting  or 
initial  phase  of  16.7  ug/As. 

Erosion  during  steady  state  operation.  During  the  long  duration 

tests,  lasting  about  240  to  60  minutes  at  steady  state  currents 

netween  1600  and  about  4000  A,  respectively,  i.e.  below  onset 
conditions  (since  steady  state  M PD  thrusters  must  be  run  below 
onset  to  avoid  destruction),  the  measured  total  mass  loss  is 
about  400  mg.  The  erosion  rate--not  counting  the  losses  of  the 
starting  phase - -amounts  to  about  0.03  ug/As  ,  which  is  three 
orders  of  magnitude  below  that  of  the  ignition  phase.  The  meas¬ 
ured  erosion  rates,  j iven  by  the  cathode  mass  loss  Am  divided  by 

the  accumulated  e  1  c  ic  charge  0  in  eg/C,  for  the  continuously 
running  steuiy  st.it**  ire  are  plotted  as  a  function  of  arc  power, 
arc  current  and  a  ;  u  t  met  ion  of  the  correlation  factor  I  /m  in 
Figs.  6ba,b,c,  respect ively,  where  m  is  the  propellant  mass 
flow.  In  all  -:■!**■  case.;  the  erosion  rate  may  be  considered  a 
e.ui.  ■:  in'  and  i  .  .  it**;  .*n  ietit  of  p,  [  and  I'  m,  i.e,  one  may  con- 

,-i  i*  i  *  •  *  lot  »  *.*•  me.::  lose  pet  unit  time,  m  _  -  Am,  At  ,  is  proper- 

V 

*  i  i  >i:  < !  *  »  •.*•  -  u  r  .  •  ■ :  i '  .  s  i  tie**  o  l  A*  ,  one  obtains  a  mass  loss 

:  it*'  ■ !  i  :■  ■  i  ’ 

.  ,r  I  ..a  ■;  ( 4  ) 


ijuas  l- s tead^  MPO  Thruster 

Measurements  on  a  pulsed  thruster  with  the  same  electrode  mate¬ 
rial  and  geometry'  show  an  average  erosion  rate  ot  about  14  ( ♦ 
3.5)  .»g  C  with  argon  as  propellant  and  16  (  t  3.5)  eg/C  with  ni¬ 
trogen  as  propellant  (see  Fig.  61).  These  latter  results  were 
obtained  by  weighing  the  cathode  before  and  after  about  420 
snots  at  the  same  operation  conditions.  The  pulsed  thruster  was 
energized  by  a  pul s e- form l ng  network  with  electrolytic  capaci¬ 
tors'0  which  produces  rectangular  current  pulses  of  up  to 
IiMux  51  8000  A  lasting  about  2  ms.  It  is  interesting  to  note 
that  tne  eros''>u  rate  measured  for  the  pulsed  thruster  is  about 
the  same  us  that,  for  the  initial  phase  of  the  continuous  thrust¬ 
er,  but  also  is  a  factor  3  to  5  greater  than  the  erosion  meas¬ 
ured  in  Japan  with  quasi-steady  thrusters .  1  3,  J u 

In  tne  following  section,  an  explanation  of  a)  the  spotty  "cold 
cathode  attachment"  with  a  fairly  high  erosion  rate  and  b)  the 
more  diffuse  "hot  cathode  attachment"  of  much  lower  erosion  will 
be  presented. 


i 


4 . 3  Theor 


Based  on  the  experimental  findings,  one  distinguishes  between 
(a)  a  spotty,  more  or  less  highly  erosive  and  (b)  a  bright  glow¬ 
ing  diffuse,  low  erosive  cathode  attachment  within  an  MPD-arc 
thruster.  First,  let  us  consider  the  spotty  arc  cathode. 


Cold  Spotty  Arc  Cathode 


The  starting  phase  of  an  MPD-arc  may  be  described  by  the  follow¬ 
ing:  at  any  favored  spot  (microprotusion  and/or  dielectric  im¬ 
purity)  on  the  cathode  surface,  the  electric  field  becomes 
strong  enough  to  start  a  breakdown.  This  breakdown  causes  a  lo¬ 
cal  overheating  of  the  emitting,  microscopic  spot  site,  thus 
leading  to  evaporation  and  a  plasma  cloud  which  expands  from  the 
spot  site  on  the  cathode  surface  into  the  gaseous  half  space. 
After  this  first  ignition  phase,  which  according  to  Guile  and 
Jiittner may  last  about  10~y  s,  a  small  crater  with  a  diameter 
of  a  few  microns  is  formed  (see  Fig.  62).  A  few  typical,  but 
more  developed  craters  left  on  the  surface  of  a  thoriated  tung¬ 
sten  cathode  are  shown  in  the  photomicrograph  of  Fig.  63. 


I’lasma  Vapor  Jet 


Fig.  b  2 
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is  much  smaller  than  the  crater  radius  rj.  Indeed,  within  the 
crater  the  plasma  pressure  and  hence  the  partial  pressure  of  the 
electron  gas  can  be  considered  as  much  higher  than  outside  ; 
therefore,  Xq  and  the  thickness  of  the  plasma  layer  above  the 
inner  crater  surface  is  thin  . .(  <Xr  -  )  while  outside  ’of f  ’the  crater 
the  plasma  layer  is  relatively  thick.  As  a  consequence  of  this 
plasma  boundary  configuration,  the  electric  field  strength  Ey  on 
the  concave  crater  surface  is  much  larger  than  outside  and 
therefore  strong  enough  to  maintain  a  high  current  density,  due 
to  field  emission  or  thermal  field  emission  on  the  inner  crater 
surface.  Tne  following  quantitative  discussion  shall  corroborate 
this  explanation. 


Fig  .  *)  4  pint  rent  density  due  to  thermal  field  emission  for  an 

electrode  work  function  t  =  4.5  eV  (tungsten). 
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•  .>•  electron  gas  of  p  =  4.6*10' 

p:  ensure  tungsten  vapor  plasma . 

■  i.;r  tie  taken  larger  than  1  0,000  K, 
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anode  :  see  Fig.  65a). 
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The  cuaeut  density  within  this  channel  is,  as  we  have  seen, 
quite  nigh  so  that  electromagnetic  forces  cannot  be  neglected. 
As  it  has  been  shown  e 1 sewher e 4 J ,  4 4  such  a  current-carrying 
plasma  channel  with  an  axial  flow  is  only  straight  and  stable  if 
the  pressure  difference  between  the  plasma  and  vapor  pressure 
pv  at  the  inner  crater  surface  and  the  ambient  pressure  p,*, 
is  larger  than  (about)  three  times  the  average  magnetic  pres¬ 
sure.  A  more  detailed  investigation  shows  that  this  factor  is 
about  3.0  for  an  almost  flat  crater  (Fig.  66);  for  a  hemispheri¬ 
cal  crater  (Fig.  65a)  it  is  about  3.3.  The  mathematical  criteri¬ 
on  for  a  stable  straight  discharge  channel  emanating  from  an  ac¬ 
tive  cathode  spot  follows  therefore  by 


p  -  p  >  k  7T—  I  a 
rv  8  »i  s  J  s 


where  3  »  k  <  3.3,  uu  =  4  n  - 1 0  ~  Vs/Am  (the  magnetic  permeabili¬ 
ty),  I3  is  the  spot  current  and  js  is  the  average  current 
density  across  the  spot  orifice.  Since  for  any  arc  spot  the  am- 
oient  pressure  p®  can  be  assumed  small  compared  to  the  spot 
pressure  pv  above  the  melt,  the  value  of  the  product  spot  cur¬ 
rent  times  current  density,  at  which  a  spot  discharge  goes  from 
stable  to  unstable,  can  be  directly  related  to  py  . 

In  Fig.  67,  the  vapor  pressure  curves  or  the  spot  conditions 
pv  and  Ts  for  stable  and  unstable  behavior  on  pure  metal 
cathodes  like  copper  molybdenum  and  tungsten  are  illustrated. 
Below  these  curves,  a  spot  discharge  is  stable;  above  these 
curves  it  is  unstable  and  cannot  exist. 


It  the  jbovf  criterion  (7)  is  not  fulfilled,  any  small  disturb¬ 
ance  leading  to  a  minute  curvature  causes  the  current-carrying 
plasma  channel  to  kink,  i.e.  the  channel  axis  bends  more  and 
more  and  eventually  the  channel  will  come  in  contact  with  the 
cathode  surface  at  a  nearby  site  off  the  original  spot  (Fig. 
65d).  The  new  site  wiLl  be  heated  up  and  a  new  current-carrying 
spot  will  be  initiated.  Again,  for  this  new  spot  the  stability 
requirement  (7)  must  be  fulfilled  in  order  to  maintain  a  stable 
and  active  discharge. 


rv.it; ,  by  inct.e  i  ;  i  ,  i  1  :n>  arc  out  reul  1  higher  and  higher,  one  has 
to  .*x;>~o"  h,  it  m . r  and  more  mierospots  are  created ,  which  now 
spread  ov-r  •  c  i*  h  «•  !-■  mi  f  ice  and,  or  cluster  together  as  larg¬ 
er  spot,  s .  TVj  -c !  a-'.  *  i :  :  is  got:  spots  can  now  have  features 

similar  to  t  v-  :n  i  or-  '.-.pots .  Again,  a  plasma  and  vapor  jet  will 
propagate  tun?  i  .-l  a  ■; ; )  spot  (Fig.  66),  and  for  the  main  cur¬ 
rent-cart  y  mg  chunite  1  wnich  emanates  from  th  is  larger  spot  site, 
the  same  ; t a n i 1 i t  y  criteria  \ 7 )  as  derived  for  the  microspot  is 
now  applicable,  oi.e  taerefore  may  state  that  the  requirement  (7) 
is  valid  very  genera  l  l y  tor  any  stable  type  of  cathode  spot.  For 
the  cl  j stored  or  1  if.K-f  spots,  however,  their  total  spot  cur¬ 
rent,  ls ,  will  h»  lirget,  while  the  average  current  density, 
j  s ,  becomes  sm a  1  I »t  •!  an  in  the  case  ot  a  microspot. 


The  vapor  pr- 

•  s  ; a r ••  a ; . 

v I*  t.  no 

;:i  e  1 1  pv 

in 

requirement  (7)  can 

now  be  relate 

■d  t  _>  "v 

«'V  jpo  r.  a  t-. 

ion  rate 

and 

to  the  average  sur- 

f  a  ce  :.  em  pe  t:  j 

...  ; 

\:ie  sijO 

t .  Hence 

,  th 

e  stability  criteria 

(  7  j  c  a  :  i  a  ;  »  : 

•  *  ’  A  !  L  t  *  *  ?  1 

ir.  f  1 1  e 

:  a  r  :r. 

ri 


(8  ) 


,i!  >-r  x  a  :  ,  ^v, iiiot  ated  per  second,  of 

it  a  ■  ■  c  we  lg at  of  the  vapor  ,  R  is 

i  I,  '< ,  i  a  geometric  factor 

i  ,.  d-qw-nd  in g  an  the  presence  of  a 
a.  :g  .  n  *  ■  i  spot  profile.  Since 
c  ,  1  i!  cei  t  han  the  ambient 
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one.  A  spot  discharge  therefore  always  requires  a  minimum  evap¬ 
oration  rate  mVs,min'  which  is  proportional  to  the  square  root 
of  the  molecular  or  atomic  weight  M  of  the  evaporated  cathode 
material  and  inversely  proportional  to  the  square  root  of  the 
melted  spot  surface  temperature,  Ts .  If  the  cathode  attachment 
consists  of  n  individual  spots  identified  by  the  indices 
1 , 2 , . . . , i , . . .  ,n  ,  the  overall  minimum  evaporation  rate  follows  by 
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and  since  the  total  current  I  is  given  by  the  sum  over  all  spot 
currents  Is,  the  overall  minimum  erosion  rate  due  to  evapora¬ 
tion,  ev  min»  1S  by  definition 


[kg/C  j. 


(  11  ) 


It  is  of  interest  to  note  that  according  to  eq.  11  the  spots 

with  relatively  high  individual  current  rates  (Isi  large), 
i.e.  the  clustered  spots,  contribute  the  major  amount  of  ero¬ 
sion,  and  that  the  erosion  rate  is  about  proportional  to  the 
spot  current  of  those  spots.  If  the  cathode  attachment  consists 
only  of  many  identical  microspots  with  the  same  average  spot 

current,  Isi  =  Is  =  const,  and  the  same  average  surface  tem¬ 
perature  Ts  =  Ts  =  const,  the  overall  erosion  rate  becomes 

independent  of  ttie  number  of  spots  and  reduces  to 
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The  smaller  therefore  the  spot  current,  the  smaller  the  erosion 
rate.  This  does  not  mean  that  the  overall  arc  current  must  be 
small,  but  rather  that  the  number  of  spots  are  large. 
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Similar  runiii^  ion;;  must  exist  now  on  the  tungsten  cathodes  of 
pulsed  thrusters,  since  during  the  short  discharge  the  cathode 
does  not  neat  up  to  thermionic  temperatures.  Ag .  .n,  spots  are 
present,  and  since  the  measured  erosion  rate  is  of  about  the 
same  amount  (14  >*g  C)  and  the  spot  temperature,  Ts ,  must  be 
taken  about  the  same  (since  the  dissipation  mechanisms  within 
spots  of  the  same  erosion  rate  are  the  same),  the  spot  current 
amounts  again  to  about  135  A  or  to  27  A  with  a  droplet  loss  fac¬ 
tor  of  5.  At  tne  current  maximum  of  8000  A,  the  number  of  spots 
are  tnerefore  between  60  and  300  .  Since  the  conditions  in  all 
these  spots  vary  slightly,  they  are  not  stable  but  should  rather 
jump  according  to  tne  following  explanation.  If  for  instance  the 
plasma  pressure  within  a  spot  increases  slightly,  the  Debye 
length  and  therefore  tne  electric  resistance  of  the  spot  de¬ 
creases  .  According  to  Kirchnoff's  distribution  law,  the  current, 
will  be  concentrated  through  that,  circuit  of  lowest  resistance-- 
m  our  case  through  tne  spot  of  lowest  resistance.  That  means, 
however,  tnat  for  tens  spot  the  right  side  of  tne  stability  re¬ 
quirement  (7),  (8)  or  (9)  increases  faster  than  the  left  side. 
Tne  consequence  will  be  that  the  discharge  channel  above  the 
spot  will  become  unstable  and,  hence,  causes  a  splitting  or  a 
jump  to  a  new  spot  site.  The  spots  are  therefore  distributed  and 
may  jump  over  the  cathode  surface,  but  according  to  inhomogene i- 
"i^s  ot  tne  catnode  surface  and  of  tne  surface  coding,  as  well 
as  the  entire  arc  snape,  they  will  prefer  to  appear  on  certain 
favorable,  still  nard  to  predict,  areas. 
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The  current  density  on  a  therm  ion ica 1 1 y  emitting  cathode  surface 
is  governed  by  the  Richardson  equation 


j  =  A0  T“J  e  kT  [ A/m  1  j  (  15  ) 

where  A  j  is  the  Richardson  constant,  t>  is  the  electron  work 
function  and  k  is  the  Boltzmann  constant.  Hence,  with  increasing 
surface  temperature,  the  current  density  increases  rapidly. 
Assuming  now  that  within  the  diffuse  attachment  area,  due  to 
small  fluctuations,  the  temperature  on  one  spot  would  increase 
only  slightly  above  the  proper  value,  the  current  density  would 
increase  in  this  spot,  and  along  with  the  current  density,  due 
to  ohmic  heating,  the  temperature  would  also  increase,  and  so 
on.  According  to  Kirchhoff's  distribution  law,  the  current  would 
be  more  and  more  concentrated  until  field  emission  sets  in  and 
self-magnetic  field  effects  play  a  role.  In  that  case,  however, 
the  above  spot  theory  becomes  applicable. 
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channel  to  touch  the  cathode  in  the  vicinity  of  the  spot  and 
eventually  leads  to  a  new  spot  site.  This  effect  causes  the  spot 
to  split  and/or  to  move  d  iscont i nuousl y  over  the  surface.  These 
spots  can  now  cluster  together  to  form  more  stationary  larger 
ones,  which  leads  to  an  enhanced  erosion,  it  is  shown  that  the 
smallest  possible  overall  erosion  rate  is  directly  proportional 
to  the  spot  current  and  that  it  is  not  necessarily  a  function  of 
the  arc  current.  Therefore,  in  order  to  avoid  a  high  erosion 
rate,  the  number  of  spots  should  be  high,  while  the  average  spot 
current  should  be  small.  The  relatively  high  cathode  erosion 
rates  measured  with  the  pulsed  and  during  the  starting  phase  of 
the  continuous  MPD-arc  thruster  are  explained  by  spots  which 
carry  a  current  of  about  30  to  160  A. 

The  processes  occuring  on  a  glowing  cathode  surface  of  a  steady 
state  continuous  MPD-arc  device  are  discussed  and,  based  on  data 
in  the  literature,  the  expected  erosion  rate  determined.  Since 
such  a  seemingly  diffuse  thermionic  attachment  is  obviously  un¬ 
stable,  it  is  suggested  that  this  so-called  thermionically  emit¬ 
ting  cathode  also  consists  of  many  tiny  hot  spots  in  which  field 
emission  is  predominant. 


5.  Numerical  Modeling  of  tne  Flow  Discharge  in  ,MPD  Thrusters 

5 . 1  Int rod uc t  ion 

Tne  increase  in  onboard  electrical  power  in  the  near  future  is 
again  raising  interest  in  self-field  MPD  thrusters.  Despite 
t.neir  simplicity  m  design  and  power  conditioning,  they  are  han¬ 
dicapped  even  today  :>y  the  shortcomings  ot  low  efficiency  and  a 
per i or, nance  limit  Known  as  “onset"  which  restricts  the  flow  rate 
at  a  given  current  value.  Therefore,  since  the  beginning  of  the 
to  s  ig  n  and  investigation  of  self-field  MPD  thrusters  at  the  end 
of  the  sixties,  e  t  fort  has  been  made  to  unterstand  the  physics 
ot  tne  MPD  discharge  and  the  development  of  appropriate  codes, 
hut  the  MPD  problems  cannot  be  solved  easily,  because  in  con¬ 
trast  to  tne  hardware  simplicity  of  these  devices,  the  physics 
involved  are  of  extreme  complexity.  To  mention  only  a  few:  real 
gas  effects,  i.e.  cold  propellant  is  heated  up  in  an  arc  and  the 
gas  has  to  be  dissociated  and  ionized;  the  problems  related  with 
tne  electrodes  are  only  partly  known  or  solved;  the  plasma  is  in 
a  thermal  non-equilibrium;  because  of  the  low  densities  and  high 
temperatures  the  Reynold's  numbers  are  low  and  the  friction  can¬ 
not  be  neglected;  the  magnetic  Reynold's  number  is  too  high  to 
separate  flow  and  discharge;  etc.  This  list  is  by  no  means  a 
complete  one.  Added  to  these  physical  complexities  are  numerical 
difficulties:  the  necessary  coupling  of  different,  partial  dif¬ 

ferential  equation  systems,  namely  elliptical  and  hyperbolic 
ones;  non- 1  inear  i  t  ies ;  steep  gradients;  and  so  on.  And  let  us 
not  forget  difficult  geometries.  Therefore,  due  to  the  great 
complexity  of  the  physical  conditions,  for  aii  codes  established 
to  solve  at  least  parts  of  these  problems,  many  and  great  sim- 
pl  i  t:  icat.  ions  must  be  made.  The  resulting  codes  include  for  exam¬ 
ple  onc-d  imens  ional  ,  '*  •  4  '  quas  i- t  wo-d  imens  iona  1  and  simplified 
two-dimensional  ones.1  , 

In  *■  h ;  s  chapter,  further  approaches  are  presented :  first,  .he 

r  projon  of  Hiigel  1  J  has  been  modified  and  used  for  different 
m'i  1  *■>•;  for  nozzle-type  thruster..;.  Second ,  a  partly  t  wo- 
.  i  1  ,  more  refined  code  has,  been  developed .  The  s  nr.pl  i- 
'  -  i  n  .ional  electromagnet  ie  calculation  fol  lows  from 


Ohm's  law,  similar  t..j  Ret.  49,  and  the  flow  is  modeled  by  a 
one-dimensional,  non-equi librium  expansion.  The  equations  are 
solved  only  for  tne  steady  phase.  The  third  approach  uses  com¬ 
parable  electromagnet  ic  equations,  but  models  the  flow  actually 
two-d imens tonal  1  y  witn  a  modified  MacCormack  code.  All  three 
approaches  neglect  el*-  '-rode  effects  and  friction. 

5.2  Ad  apt  at  ion  of  tne  Semi-Fmp  i  r  ica  1  Hiigel  Calculation  Method 
to  Include  Different  Geoinetries 

In  Ref.  10  Huge!  developed  a  magnetof 1 uiddynam ic  model  which  he 
used  to  explain  the  nature  of  the  discharge  instabilities  set 
ting  in  at  the  so-called  critical  value  of  I  V ni .  In  this  model 
(Fig.  63),  tne  similarity  between  the  measured  flow  field  with  a 
spatial  source  flow  field  is  used  to  establish  a  computational 
code  to  resolve  the  midstream  flow  in  a  quas i- two-d imens iona 1 
manner.  As  a  further  assumption,  the  magnetic  Reynolds  number  is 
taken  as  low,  i.e.  Rm  <  1.  In  such  a  case,  it  is  possible  to  de¬ 
termine  the  electromagnetic  field  separately  from  the  flow 
field,  whereas  its  effect  upon  the  latter  is  still  explicitly 
taken  into  account  via  force  and  heating  terms  in  the  momentum 
and  energy  equal-  ions .  As  a  first  step,  then,  the  distribution  of 


8 1  - 


current  density  and  magnetic  field  is  calculated  i >y  using  poten- 
ttal  theory  and  Maxwell's  equations  and  by  introducing  a  stream 
function  -K  The  electrode  geometry  and  tin*  current  distributions 
on  the  electrodes,  needed  as  boundary  conditions  for  the  solu¬ 
tion  of  the  partial  differential  equation  for  v,  are  taken  from 
the  experiments.  Good  agreement  between  calculated  and  measured 
current  distributions  justifies  this  procedure. 

Hence,  from  the  usual  set  of  basic  equations,  i  .  e .  state,  conti¬ 
nuity,  momentum,  energy,  and  Ohm's  law,  a  system  of  partial  dif¬ 
ferential  equations  for  the  flow  properties  is  obtained.  Jo u  1  e 
heating  and  electromagnetic  forces  are  included,  whereas  heat 
conduction,  radiation,  and  viscous  effects  are  neglected.  In 
particular,  the  results  of  that  model  indicated  a  reduction  of 
the  pressure  at  the  anode  wall  witn  increasing  I  *7  m .  The  value 
at  which  the  pressure  approached  zero  agreed  well  with  the  ex¬ 
perimentally  observed  value.  Although  a  wide  variety  of  parame¬ 
ters  had  been  investigated,  both  experimentally  and  theoreti¬ 
cally,  the  geometry  of  tne  thruster  had  never  been  altered. 
Hence,  it  appeared  necessary  to  also  cover  this  aspect  in  the 
research  work. 

Along  this  line,  then,  a  study  was  made  to  prove  whether  this 
magnetof 1 uiddynam  ic  model  could  be  used  as  an  appropriate  tool 
for  parametric  and  predictory  calculations.  '  *  Thereby,  the  main 
goal  was  to  investigate  the  effect  of  the  opening  angle  of  the 
thruster's  cone  nozzle  upon  the  critical  conditions  and  the 
achievable  thrust. 


The  original  model  was  changed  to  omit  some  s imp!  i  f  icat ions  ap¬ 
plied  in  the  mathematic  treatment  of  the  cone  geometry.  In  par¬ 
ticular,  when  calculating  volume  elements,  stream  functions  and 
electromagnetic  force;  smal  l-ang  l**-approx  imat  10ns  when  using  go- 
niometric  functions  were  replaced  by  the  exact  expressions.  In 
Fig.  69,  computed  value;;  of  tne  pressure  at  the  anode  are  plot¬ 
ter!  as  a  function  of  tne  cor re  1  at  ion  parameter  1  m  for  differ¬ 
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"critical"  condition  Clearly,  tnis  appears  as  a  consequence  of 
tne  fact,  that  with  increasing  1 ‘ ,  m  the  self-magnetic  forces 
lead  to  a  steepening  pressure  profile  over  the  cross  section. 
Note  that  (lVm)crit  is  the  same  for  all  values  of  in. 

The  effect  of  the  opening  angle  upon  the  critical  value 
(lVrn)Crit  i-s  summarized  in  Fig.  70.  It  appears  that  some 
asymptotic  value  of  about  4  *  1 0 1 °  As "/kg  is  reached.  From  these 
results,  it  should  be  advantageous  to  use  thrusters  with  large 
open  ing  ang les  . 

However,  of  much  more  interest  is  the  thrust  that  can  be 
achieved  at  the  corresponding  conditions.  Or  more  prec i se 1 y , can 
one  expect  an  increase  in  specific  impulse?  In  Fig.  71,  the  cal¬ 
culated  thrust  at  constant  mass  flow  rate  is  snown  as  a  function 
of  current  and  opening  angle.  The  dashed  curve  represents  the 
limits  imposed  by  reaching  the  critical  values  of  ll/m.  The 
aasned  curve  "3b°eXp"  show  the  experimental  results  for  the 
case  of  an  opening  angle  of  36°  (Stuttgart  thruster)  which  coin¬ 
cide  quite  well  with  the  calculated  values.  There  is  strong  evi¬ 
dence  that  the  specific  impulse  rather  decreases  than  rises  with 
increasing  angle. 


iz  «  so  »c  '0 

a  - — > 


K  i  |  .  'U  on:;**'  >  •  r  i  t  ♦*  r  1 « >n  .in  .  i  f  1 1  rnM  i  on  * » 1  t  :  1  •  ■  <  >p<  ■  n  1  nj  .  i  nj  1  . 


v\  W**> 


"j .  3  Steady  -  State  _MPD  T  h_r  us  t  er  1  n  v  es  t  lg  at  ions 

In  order  to  predict  the  overaLl  performance  of  continuously  run¬ 
ning,  self-field  MPD  thrusters,  a  sem  i- 1  wo-d  linens  iona  1  model 
calculation  and  computer  program  has  been  developed.’1,^11  In  a 
first  step  trus  program  has  been  used  to  investigate  the  current 
contour  lines  and  the  thrust  within  the  nozzle-type  thruster  DT2 
(Fig.  3). 

Current  Contour  Lines 

In  this  approach,  the  current  contour  lines  are  calculated  in  r- 
and  z-  dependence  while  the  flow  field  is  taken  as  a  one-dimen- 
sional,  frictionless  expansion  flow,  assuming  isothermal  behav¬ 
ior  for  the  electron  component  and  isentropic  expansion  tor  the 
heavy  particles  like  ions  and  neutrals.  This  assumption  is  jus¬ 
tified,  since  primarily  the  electron  component  is  heated  by 
Ohmic  energy  input,  and  since  within  the  nozzle  expansion  flow 
the  energy  exchange  between  the  lightweight  electrons  and  the 
heavy  particles  like  ions  and  neutrals  is  fairly  weak.  There¬ 
fore,  reaction  collisions  (ionization  and  recombination  process¬ 
es,  etc.)  are  neglected.  Hence,  the  model  approach  distinguishes 
between  electron  and  heavy  particle  temperatures  and  assumes  a 
frozen  flow  from  the  nozzle  throat  on  downstream. 

Within  the  nozzle  throat,  one  assumes  a  "fully  developed"  arc 
flow  regime,  i.e.  the  gas  and  plasma  velocity  and  the  current 
density  vector  have  only  an  axial  component  and  the  arc  column 
is  defined  by  only  radial  heat  transfer  (and  possibly  by  radial 
radiation)  losses.  Moreover,  from  the  nozzle  throat  on  upstream, 
one  no  longer  distinguishes  between  electron  and  heavy  particle 
temperatures,  since  within  this  higher  pressure  regime  a  tight 
coupling  between  all  components  is  more  likely  to  be  established 
than  in  the  low  pressure  downstream  nozzle  regime.  In  spite  of 
the  neglect  of  the  radial  dependency  of  the  flow  and  species 
temperatures,  however,  this  model  approach  turns  out  to  be  an 
effective  tool  for  predicting  the  current  contour  lines  within 
and  the  thrust  and  specific  impulse  of  a  nozzle-type  self-field 
MPD  th ruster  . 


Assum 


•  steady  state  conditions 

•  rotational  symmetry 

•  no  azimuthal  current 

•  singly  ionized,  quasi-neutral  plasma 

•  electric  current  density  normal  to  electrode  surfaces 

•  from  cathode  to  nozzle  throat:  thermal  equilibrium 

•  within  the  nozzle:  frozen  flow,  ambipolar  expansion  flow 

without  friction,  electrons:  isothermal,  ions:  adiabatic 

•  outside  the  nozzle:  hyperbolic  expansion  flow,  electrons 

and  ions  adiabatic 


Basic  Equations 
Maxwell's  eqs .  : 


Ohm  s  haw : 
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State: 
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Continuity  Equation: 


;>vA  =  o  A  ^  v ,  =  m 
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Expansion  Relation: 
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Nozzle  Throat  Conditions: 
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Tt  =  Te  =  f  •  (|-)2/5  (25) 

where  f  is  a  slowly  varying  function  depending  on  current  dis¬ 
tribution  . 

The  last  three  equations  are  derived  in  detail  under  the  above 
assumptions  in  Ref.  4. 

Rewriting  Ohm's  Law  by  means  of  Maxwell's  equations,  one  obtains 
a  vector  equation  for  the  magnetic  induction  vector  E?  in  the 
form 


_5B 


N 


rot 


1  * , 

1  -  rot  B  I 
o 


rot  | v  x  B j 


+  rot  P  (rot  5)  x  E?  1  =  0 

a  u 


(26  ) 


with 


.)  t  _  1 

B  o  en 


(27  ) 


For  a  stream  function 

-  rB0 

with  respect  to  the  rotational  symmetry  and  the  zero  azimuthal 
current,  the  elliptic,  partial  differential  equation  of  2nd  or¬ 
der  follows  from  eq .  (26)  in  the  form 
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The  function  f(r,z)  =  const,  represents  now  a  current  contour 
line,  since  D  =  8.-.  "  I  (  r )  r  and  1  (  r )  is  the  electric  current 
carried  through  a  cross  sectional  area  of  ir  .  The  proper  bound¬ 
ary  conditions  for  ;  follow  from  the  geometry  of  the  thruster 
wa 11s  and  e 1 e c t  rode s  . 


A  solution  of  tne  equation  system  (16)  -  (28)  requires  as  input 

data: 

•  Geometry  of  Thruster  A(z) 

•  Mass  Flow  Rate  m 

•  Electric  Current  I 


and  yields  tne  following  results: 

•  Current  Contour  Lines 

•  Thrust  f 

•  Specific  Impulse  iSp 

Step  1:  By  setting  t  no  integral 
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in  the  expansion  relation  eq .  (23)  to  zero,  one  obtains 

together  with  ttie  nozzle  throat  conditions  (24)  and 
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which  means 

v  r  *  0  (Fig.  72  )  . 
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knowledge  of  t:ie  flow  ana  plasma  properties 
(m  first  approximation),  ',  .■  and  v  are  also  known,  as 
well  as  their  derivatives.  With  these  quantities  one 
now  solves  the  differential  equation  28  by  means  of  a 
finite  differential  metnod  { Gauss-Se ide 1 )  for  the  cur¬ 
rent  contour  lines  in  first  approximation.  In  Fig.  73 
trie  current  contour  lines  are  shown  for  different  cur¬ 
rents  between  1.3  and  4  kA.  According  to  Maxwell's 
equation,  one  also  obtains  from  r  the  current  density 
ana  tne  j  x  R  force  configuration  within  the  discharge. 

Step  2:  One  solves  tne  expansion  relation  ( eq .  23)  numerically 
uy  tukini  ( j  x  3)  within  the  integral  from  Step  1. 
With  that  tne  velocity  v (  z )  and  the  quantities  p (  z )  , 

.  i  zj  ,  n,y  z  are  calculated  in  second  order  approx  lin  a- 
tion  and  from  tnat  the  new  quantities  for  o,  ?,  and  v 
with  t.neir  derivatives. 

With  these  new  quantities,  one  now  solves  the  differen¬ 
tial  equation  (28)  for  the  current  contour  lines  by 
means  of  tne  same  numerical  method  as  in  Step  1. 
Repeating  Step  2  yields  an  iterative  solution  of  tne 
equation  system  (16)  -  (25). 
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Fig.  72  Scheme  of  the  density  model.  Inside  the  nozzle:  one- 

dimensional.  Outside  the  nozzle:  elliptical. 
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Tne  calculation  shows  that  after  two  steps  the  current 
contour  lines  with  I  -  4  kA  do  not  change  anymore  and 
even  a  first  step  approach  -jives  -jood  results.  For 
lower  current  rates  the  deviation  is  even  smaller. 

A  comparison  between  calculated  and  measured  current  pattern 
(Fig.  74)  shows  a  fairly  good  agreement.  The  fact  that  in  the 
measured  case  the  anode  attachment  is  more  concentrated  down¬ 


stream  of  the  nozzle-  end  than  in  the  calculated  case  may  be  at¬ 
tributed  to  the  theoretical  model  which  neglects  the  boundary 
layer  along  tne  nozzle  and  anode  walls. 
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Fig.  'M  a ;  Calculated  and  b)  measured  current  contour  lines  of 

tne  .MFC  discharge _ in  the  nozzle-type  DT2-IRS  thruster  for  I  = 

■  ,  )  KA  and  m  -  0.8  j_  argon. 
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According  to  tiio  known  j  x  B  force  configuration  and  the  flow 
conditions,  one  can  now  calculate  the  thrust. 

The  tnrust  of  a  self- field  MPD  accelerator  is  composed  of  two 
effects;  first  the  sum  of  all  gas-dynamic  surface  forces  on  the 
inner  plenum  chamber  and  on  the  inner  nozzle  wall,  and  second 
the  sum  of  all  volume  effects  given  by  the  magnetic  field  forces 
which  act  between  the  current-carrying  plasma  and  the  entire 
power  circuit,  including  the  coaxial  electrodes.  Hence  it  is 
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where  ■  cvv  +•  P-  is  the  gas  dynamic  force  per  unit  area;  Ag  rep¬ 
resents  tne  surface  of  all  internal  walls  like  plenum  chamber 
with  the  gas  inlet  cross  sections,  the  cathode  surface  and  the 
inner  solid  wall  of  the  nozzle;  V  is  the  current-carrying  vol¬ 
ume.  The  second  term  is  the  sum  over  all  {  x  B  forces  which  act 
between  the  current-carrying  solid  leads  of  the  electric  power 
circuit  and  the  current-carrying  plasma;  in  the  second  equation 
p  is  the  chamber  pressure,  At  the  nozzle  throat  area  and 
cxn  tne  thrust  coefficient  matching  the  gas  dynamic  effects. 
It  is  by  definition 

c  =  -r — r—  j  I  dvv  +  P  f  *dA  (  10  ' 

I  P  ,  A  . 
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Tne  calculated  thrust,  based  on  the  current  density  and  mag m  *  . 
field  determined  wits  tne  procedure  as  described  in  tne  :  > 
"Current  Contour  linos" ,  i  now  compared  with  the  ex;-  ■  .  ■ 
results  in  Pin.  7:>. 
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Fig.  7  5  .Comparison  between  calculated  and  measured  thrust  of 
the  DT2-IRS  thruster.  Tm/tnax  is  the  maximum  possible  electro¬ 
magnetic  thrust. 
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5 . 4  Time  Dependent,  Two-Dimensional  MPD  Thruster  Investigations 

In  addition  to  the  steady  state  MPD  theory,  a  real  two-dimen¬ 
sional  numerical  code  has  been  developed. b4  This  code  allows  the 
calculation  of  the  time  dependent  plasma  properties  during  the 
starting  phase  of  steady  state  MPD  thrusters,  as  well  as  for 
pulsed  mode  thrusters.  The  numerical  solution  is  done  in  a  man¬ 
ner  similar  to  that  of  Jacobi bb  for  laser  applications.  In  a 
first  approach,  this  MPD  code  was  applied  to  the  cylindrical  MPD 
thruster  ZT1  being  developed  at  this  institute.  Fig.  76  shows 
the  drawing  of  ZT1,  and  Fig.  77  shows  the  corresponding  simpli¬ 
fied  scheme  for  the  calculation. 


Fig.  76  Cylindrical  continuous  MPD  thruster  ZT1-IRS. 


Fig.  77  Simplified  scheme  used  for  the  calculation, 


The  geometrical  configuration  on  which  the  computation  is  based 
was  chosen  as  close  as  possible  to  the  experimental  device.  The 
radii  of  cathode  and  anode  are  taken  according  to  the  dimensions 
of  the  real  device;  however,  the  length  of  the  cathode  is  as¬ 
sumed  to  be  of  the  same  length  as  the  anode  and  it  matches  the 
length  of  the  three  anode  segments  (see  Fig.  76).  The  model  is 
assumed  to  be  an  infinite  tube  of  an  insulator  with  a  central 
insulating  rod,  only  the  electrodes  will  be  of  zero  resistance. 

For  the  time  dependent  calculation  of  the  MPD  flow,  two  inde¬ 
pendent  codes  were  correlated  to  determine  flow  field,  current 
and  magnetic  field  distribution.  These  two  codes  are  connected 
in  the  following  manner:  for  a  given  flow  field,  the  current, 
and  hence  the  magnetic  field  distribution,  was  determined,  with 
these  results  the  flow  field  equations  were  integrated.  In  the 
next  time  step,  this  new  flow  field  and  the  new  distribution  of 
the  thermal  properties  were  taken  to  calculate  the  new  electro¬ 
magnetic  field  distribution,  and  so  on. 

As  input  data  the  two  codes  require: 

•  geometry  of  the  thruster 

•  mass  flow  rate 

•  electric  current  as  a  function  of  time 

and  at  the  starting  point: 

•  temperature  field 

•  Mach  number. 

This  yields  the  following  time  dependent  results: 

•  flow 

•  temperature 

•  pressure 

•  Mach  number 

•  magnetic  field 

•  current  density  distribution  within  the  channel 


Electromagnetic  Code 


The  electromagnetic  properties  are  obtained  in  a  similar  manner 
as  from  the  steady  state  thruster  code,  using  the  following  as- 
sumpt  ions  : 

•  rotational  symmetry 

•  no  azimuthal  current 

•  electric  field  normal  to  the  electrode  surfaces. 

One  solves  Ohm's  law  with  respect  to  Maxwell's  equations;  this 
means  one  must  solve  the  elliptic  partial  differential  equation 
28. 

The  boundary  conditions  are  as  follows: 

On  the  insulator  surface  and  at  the  boundary  of  the  integration 
area,  the  stream  function  ’F  equals  zero  for  the  case  of  the 
downstream  insulator  and  equals  ^max  for  the  upstream  insula¬ 
tor.  At  the  electrodes  the  boundary  condition  follows  from  the 
requirement  that  the  electric  field  S  must  be  perpendicular  to 
the  surface.  This  means  for  a  cylindrical  channel  as  considered 
here  that  at  the  electrode  surfaces  Ez  equals  zero.  Therefore, 
the  boundary  condition  follows  to 

J —  i!  _  vr  -  -  JLL-  ^  =  o  (3i) 

an0r  dr  r  ^0r2  dz 

With  this  closed  boundary  condition,  the  computation  was  done 
with  a  Gauss-Seidel  code  and  yields  as  a  result  the  magnetic 
field  and  the  current  density  distribution  within  the  channel. 

Flow  Field  Code 


For  the  description  of  the  time-dependent,  two-dimensional, 
supersonic,  ohmically  heated  flow  the  conservation  equations  for 
mass,  impulse  and  energy  were  used  with  the  following  assump- 
t ions : 

•  friction  free  flow 

•  rotational  symmetry 

•  no  azimuthal  current  and  velocity 

•  supersonic  flow  within  the  channel 
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Tne  conservation  theorem  formulation  was  chosen  in  order  to  sat¬ 
isfy  the  Rankin-Hugon iot  shock  conditions  implicity;  hence  the 
description  of  shock  depth  and  shock  velocity  are  very  good. 
These  considerations  yield  the  following  non-linear  hyperbolic 
differential  equation  system  with  cylindrical  coordinates 
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where  the  indices  I.,  n  ,  n  indicate  the  partial  differentia- 

t  r  z 

tion  with  respect  to  time  and  to  the  r-  and  z-directions . 
Boundary  Conditions: 

upstream:  at  that  boundary  all  flow  parameters  are  assumed 

to  be  constant  and  must  be  set. 

downstream:  here  the  flow  parameters  are  extrapolated  because 

of  the  calculated  field  parameters, 
at  the  side:  no  normal  velocity;  normal  differentiations  equal 
zero . 

Initial  conditions: 

At  t  =  0  the  flow  parameters  equal  the  undisturbed 

parallel  flow; 

uniform  current  density; 

Mach  number  and  temperature  must  be  set. 

For  the  numerical  solution  a  modified  MacCormack  code  was  used; 
this  is  an  explicit  code  of  second  order  which  uses  the  shock 
capturing  method.  It  allows  a  fairly  high  energy  input  into  the 
flow  due  to  ohmic  heating  and  a  steep  rise  in  the  impulse  due  to 
j  x  ^  forces.  Because  of  instabilities  at  the  starting  phase 
caused  by  initial  and  boundary  conditions  on  the  upstream  bound¬ 
ary,  the  initial  conditions  will  be  restricted  for  the  moment  to 
values  which  result  in  cowling  numbers  smaller  than  1.5  in  the 
middle  of  the  electrodes. 


Starting  with  radially  parallel  and  homogeneous  current  contour 
lines  at  t  =  0,  the  current  is  driven  more  and  more  downstream — 
stronger  near  the  anode  than  at  the  cathode--and  bulges  at  the 
electrode  ends.  The  dotted  lines  are  those  at  20  us  after  start¬ 
ing,  and  the  dashed  and  solid  lines  are  for  160  us  and  200  us, 
respect ively . 

Fig.  80  illustrates  the  temperature  distribution  100  us  and  200 
us  after  starting.  The  strong  rise  in  temperature  near  the 
electrode  is  caused  by  ohmic  heating  due  to  the  relatively  high 
current  density  there.  The  disturbance  at  the  cathode  end  due  to 
discontinuous  boundary  conditions  decreases  with  time.  This 
fairly  marked  disturbance  at  starting  is  one  reason  for  choosing 
a  high  initial  temperature. 
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Fig.  80  Temperature  distribution  within  the  channel  at  a)  1 
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A  radial  shock  front  formed  at  the  end  of  the  electrodes  moves 


upstream  while  lowering  its  intensity;  this  is  demonstrated  in 
Fig .  81  . 


The  axial  and  radial  velocity  distribution  at  200  ns  are  plotted 
in  Fig.  82.  At  both  ends  of  the  cathode  there  is  a  steep  rise  in 
the  radial  velocity  vr,  while  the  axial  velocity  reaches  its 
maximum  at  the  beginning  of  the  cathode.  The  axial  velocity  de¬ 
creases  towards  the  cathode  end,  with  a  low  rise  from  this  point 


downstream  as  well  as  towards  the  anode. 


The  relationship  of  the  Mach  number  distribution  to  the  tempera¬ 
ture  and  velocity  distribution  at  200  ns  is  shown  in  Fig.  83. 
The  Mach  number  set  equal  to  1.2  at  the  upstream  boundary  reach¬ 
es  an  initial  maximum  at  the  beginning  of  the  cathode,  rises 
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Fig.  81  Pressure  distribution  within  the  channel  at  a)  100  ns, 


b)  200  ns. 


o.so 


\/ z  •  b)  radial  velocity  vr . 


Fig.  83  Mach  number  distribution  at  200 
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ayain  while  reaching  an  interim  maximum  at  about  1/3  of  the 
cathode  length  and  rises  once  more  behind  the  electrodes.  The 
minimum  at  the  end  of  the  cathode  is  caused  by  the  temperature 
instability  at  that  point.  It  is  remarkable  that  the  Mach  number 
is  below  1  within  a  wide  range,  especially  before  the  end  of  the 
electrodes . 

As  described  in  section  5.3,  the  thrust  is  composed  of  an  MPD 
part  and  a  gas  dynamic  part  (see  eq.  29a).  The  pure  MPD  thrust 
can  be  calculated  by  an  integration  of  the  j  x  B  forces  within 

the  entire  current-carrying  volume;  another  possibility  is  to 

1  2 

use  the  well-known  formula 


where  r./r  is  the  anode  to  cathode  ratio.  Both  have  been  done 

r\  C 

for  the  example  described  here  with  a  coincidence  of  within  1%. 


In  Fig.  84  the  pure  MPD  thrust  (dashed  line)  and  the  total 
thrust  (solid  line)  are  plotted  depending  on  time.  The  total 
calculated  thrust  at  400  ns  lies  within  a  3%  deviation  from  the 
expected  thrust. 


Fig.  84  Tnrust  versus  time.  Solid  line:  total  thrust.  Dashed 

line:  electromagnetic  thrust. 
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Parallel  to  the  MPD  investigations,  some  mission  analysis  was 
performed  to  evaluate  the  employment  limits  and  to  justify  the 
development  of  those  thrusters.  The  output  was  first  an  optimi¬ 
zation  study  in  which  the  known  results  of  Langmuir  and  Irving 
[56,  57 J  were  refined  by  better  and  more  realistic  assumptions, 
e.g.  by  taking  the  specific  power  of  the  power  supply  unit  as  a 
function  of  the  specific  impulse.  The  results  were  presented  at 
the  21st  Joint  Propulsion  Conference,  1985  (see  chapter  7)  and 
will  be  published  in  summer  1987  in  the  Journal  of  Propulsion 
and  Power. 

Second,  a  trajectory  optimization  program  [58]  was  modified  to 
allow  the  calculations  of  optimal  low  thrust  trajectories.  This 
study  is  also  supported  by  an  ESA  contract  and  is  not  yet  com¬ 
pleted  (see  Samland,  chapter  7). 
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